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ABSTRACT 

The Advanced Camera for Surveys (ACS) Fornax Cluster Survey is a Hubble Space Telescope pro- 
gram to image 43 early-type galaxies in the Fornax cluster, using the F475W and F850LP bandpasses 
of the ACS. We employ both ID and 2D techniques to characterize the properties of the stellar nuclei 
in these galaxies, defined as the central "luminosity excesses" relative to a Sersic model fitted to the 
underlying host. We find 72 ± 13% of our sample (31 galaxies) to be nucleated, with only three of the 
nuclei offset by more than 0'.'5 from their galaxy photocenter, and with the majority of nuclei having 
colors bluer than their hosts. The nuclei are observed to be larger, and brighter, than typical Fornax 
globular clusters, and to follow different structural scaling relations. A comparison of our results to 
those from the ACS Virgo Cluster Survey reveals striking similarities in the properties of the nuclei 
belonging to these different environments. We briefly review a variety of proposed formation mod- 
els and conclude that, for the low- mass galaxies in our sample, the most important mechanism for 
nucleus growth is probably infall of star clusters through dynamical friction, while for higher mass 
galaxies, gas accretion triggered by mergers, accretions and tidal torques is likely to dominate, with 
the relative importance of these two processes varying smoothly as a function of galaxy mass. Some 
intermediate-mass galaxies in our sample show a complexity in their inner structure that may be the 
signature of "hybrid nuclei" that arose through parallel formation channels. 

Subject headings: galaxies: clusters: individual (Fornax, Virgo); galaxies: elliptical and lenticular, cD; 
galaxies: nuclei; galaxies: structure 



1. INTRODUCTION 

Once viewed as relatively simple objects that formed 
in a single, "monolithic" collapse, early-type galaxies 
are now widely believed to have been assembled hierar- 
chically through repeated mergers and accretions (e.g., 
White 'fe Reeslll978tlSearle fe ZiniJ[l97l iWhite fe Frenkl 
19911: iKauffmann fe Haehneltl 120001: iCole et al.l 120001: 



Springel et ail 120051: iBower et al1l2006l ). A property of 
most luminous (e.g., M r < —22.5) early-type galaxies is 
that they appear to have formed the majority of their 
stars at high redshift (z > 1, corresponding to ages of 
t > 7-8 Gyr) and on short timescales (At < 1 Gyr) 
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/e.g.. IBower et al.lp9l lFranxj|j~9 93: Tho mas et aHll99l 
iTrager et all 120001 : iWake et al.l l2006tT These features 
may be related to feedback from active galactic nuclei 
( AGN) , which can generate jets and outflo ws that blow 
away gas a nd su ppress star for matio n fe.g.. iSilk fe Reesl 
19981: lKindl20rjl iMurrav et alJl^OCjl iFabian et al.ll2006t 
Robertson et al.l l2006( ). The general trends in the star 



formation histories of low- and intermediate-luminosity 
early-type galaxies are not as well understood, but they 
are known to show considerable div ersity and to depen d 
sensitively on environment (see, e.g.. lTolstov et al"T 2009). 

The discovery of th e Mr^-a r elation 

(jFerrarese fe Merrittl [20001 IGebhardt et~aTl [H(| 
points to a fundamental connection between the central 
black holes powering these AGN, and the dynamical 
properties of their host galaxies. There are several 
other galaxy properties that have also been found 
to sc ale with black hole mass, i ncluding luminosity 
(e.g.. iKormendy fe Richstond 119951: IFerrarese fe Merrittl 
l2000p . light concentration Ce.g.. IGraham et al.l I2001D 



global velocity dispersion (e.g., Ferrarese fe Merrittl 



120001: IGebhardt et ail 120001 : iGultekin et al.l 120091) bulge 
(e.g..lMagorrian et al.lll998l : lMarconi fe Hundl20q 



mass 



IHarine fc RbJl2 004L and t otal gravitational mass of the 
host (jBandar a et al.1 120091 ) . Thus, it has become clear 
that an understanding of the central regions of galaxies, 
including black holes and AGN, is essential if we are to 
make sense of the formation and evolution of galaxies 
themselves. 

However, the direct detection of black holes re- 
mains very cha ll engin g: see, e.g., Chapter 11 of 
IFerrarese fe Farcfl ([2005T ) for an overview of the observa- 
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tional difficulties. For kinematic measurements, a high 
central surface brightness is needed to obtain spectra 
of adequate S/N, and this requirement can pose prob- 
lems for massive early-type galaxies with shallow surface 
brightness profiles in their cores. At the distances of the 
Virgo and Fornax clusters, the small angular size of the 
black hole "sphere of influence" in most galaxies intro- 
duces a further complication. For example, at 20 Mpc, 
the distance of Fornax, a black hole in a galaxy with 
a = 200 km s _1 has a sphere of influenc e of only 0'/2 in 



radius (assuming the A4-a relation from iFerrarese et al.l 
2006b]). It is therefore not surprising that a dynami- 
cal black hole mass measurement exists for only a sin- 
gle early-type galaxy in the Fornax cluste r (FCC 213; 
iHoughton etalll200a iGebhardt et al.ll2007ft . 

On the other hand, the correlation between a 
galaxy's mass and that of its black hole was re- 
cently shown to extend down to the central nuclear 
star cl u sters found in low - mass galaxies (|Ferrarese et al.l 
I2006al iWehner fe Harris! I2006D . Other studies have 
reported similar relationships between black hole or 
nucleus mass an d the host bulge luminosity, mass , 
and Sersic index (iRossa et al.ll2006t iBalcells et alJl2007t 
IGraham fc Driverll2007t) .' These results are suggestive of 
a global relationship between galaxies and bot h types of 
central massive object (CMO: lC6te et al.ll2006l hereafter 
C06): however, it is still an open question as to whether 
black holes and nuclei form via the same mechanisms, 
or whether nuclei form first and serve as seeds for black 
hole formation. 

The hydrodynamical simulations of iLi et al.l (|2007f ) 
of a shared formation mechanism for both nuclei and 
black holes via the gravitational collapse of gas in bul- 
gelcss disks were able to reproduce a CMO and host 
mass correlation even without imposing an a priori M- 
a relation, and were o bserved to be in agreement with 
IFerrarese et al.l ()200 6a). Alternatively, IFerrarese et al.l 
(2006a) noted that nuclei could, in principle, form in 
all galaxies, but in massive galaxies they might either 
collapse or be destroyed (or otherwise altered) by bi- 
nary black holes. U sing semi-analytic mo d els, it was 
the demonstra t ed by iDevecchi fc Volonteril (|2009l ) and 
iDevecchi et al.l (|2010l ) that nuclei could form at high red- 
shifts and act as possible black hole seeds. 

If nuclei and black holes form simultaneously, then 
it is possible that momentum feedback determines 
which object will eventu ally dominate the CMO mass. 
iMcLaughlin et al.l (|2006l ) noted that the same momen- 
tum fluxthat drives out gas from black holes ([Kind 
120031 12005ft co uld also regulate the g rowth of nuclear 
star clusters. iNavakshin et al.l (|2009t ) used this find- 
ing to explain why nuclei, not black holes, appear more 
likely to form in less massive hosts. Both objects can 
form simultaneously as gas is driven to the center of a 
galaxy through an event such as a merger, but it is the 
mass of the host bulge that sets the individual forma- 
tion rates. Some evidence for such a scenario comes 
from observations of intermediate-luminosity galaxies 
(iFilippenko "fcHol [200l iGonzalez Dekado et all [200l 
iSeth et al.l l2008t IGraham fc Spitler! 120091). as we ll as a 
number of dwarfs (jBarth et al.ll2004l ; TReines et al.ll2011l ). 
that have been found to contain both a central stellar 
nucleus and a blac k hole. Indeed, u sing observations in 
the Virgo cluster, iGallo et ail (|2010l ) estimated that hy- 



brid nuclei could occur in 0.3-7% of galaxies with stellar 
masses below 10 11 M Q , and in less than 32% of hosts 
above this stellar mass. In short, the study of nuclei 
presents us with a new opportunity to deepen our un- 
derstanding of how galaxies and black holes co-evolve. 

Like black holes, nuclei pose some observational chal- 
lenges of their own. Although their existence in some 
dwarf galaxies has been known for decades, comprehen- 
sive surveys of galaxy clusters — in which the frequency 
of nucleation within complete galaxy sa mples could be 
robust ly measured — did not appear until iBinggeli et aD 
(1987) published their Virgo Cluster Catalog (VCC). 
This program observed 1277 members and 574 probable 
members of the Virgo cluster using the 2.5 m Las Cam- 
panas telescope; about 26% of all dwarf galaxies in the 
VCC sample were found to be nucleated. Shor tly there- 
after, a similar survey of the Fornax cluster by iFergusonl 
(1989) — the Fornax Cluster Catalog (FCC) — found 
nuclei in 103/249 « 41% of their dwarf galaxies. In the 
above studies, dwarf galaxies were identified primarily 
morphologically by their flat surface brightness profiles, 
although in gene ral they were found to be fainter than 
M B ^ -18 mag ()Sandage fc BinggelillT98l . 

Given the low luminosities and small sizes of most of 
these nuclei, the frequencies of nucleation estimated from 
ground-based photo graphic stu dies are certainly lower 
limits. For instance. lLotz et al.l (|2004D used WFPC2 on 
the Hubble Space Telescope (HST) to observe 69 dwarf 
elliptical galaxies in both Virgo and Fornax, finding nu- 
clei in six galaxies that were previously classified as non- 
nucleated in the VCC and FCC. Based on wide-field 
imagi ng of Virgo dwarfs from the Isaac Newton Tele- 
scope, IGrant et al.l (2005) was able to identify many faint 
nuclei that were missed in the earlier photographic sur- 
vey. In fact, imaging of late-type galaxies with HST com- 
monly revealed "nuclear clusters" that had gone unno- 
ticed in earlier studie s, with an overall frequency of nucle- 
ation of ^ 70% (e g..lCarollo et al.lll998tlMatthews et al 
1999t iBoker etaTll2004l: iWalcher et aTTl2005l; ISeth et al 



2006). 



The first study to find a comparable frequency of nucle- 
ation among early-type galaxies was carried out by IC06I 
with the Advanced Camera for Surveys (AC S) on HST: 
i.e., t he ACS Virgo Cluster Survey fACSVCS: IC6te et al.l 
2004)0 In addition to establishing a high frequency of 
nucleation for early-type galaxies (at least 66% for galax- 
ies brighter than Mb ~ —15), the high- resolution imag- 
ing made it possible to characterize the detailed proper- 
ties of the nuclei for the first time, including their lumi- 
nosity function, structural properties, color-magnitude 
relation, and nucleus-to-galaxy luminosity ratio. We 
note here that although in C06 and this work, we call 
the central excess of light rising above a galaxy's ex- 
trapolated outer surface brightness profile a "nucleus", 
these objects are not limited to being nuclear star clus- 
ters; certainly, some could be described as disks, bars, or 
other large scale structures, which have been observed 
by previous studies of early-type galaxy centers (e.g. 
Ferrarese et~aTll2006bt TBalcells et aLll2007t iMorelli et al.l 



2010). In this paper, which is part of the ACS Fornax 



12 Related papers from t he ACSVCS on the centr al structure of 
early- type galaxies include Ferrarese et al. (2006b a); Cote et al. 



2007); Glass ct al. (2011). 
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Cluster Survey (ACSFCS), we examine the properties of 
nuclei belonging to galaxies in the Fornax Cluster, which 
is located at a distance of D = 20 ± 0.3 ± 1 . 4 Mp c (sta- 
tistical + systematic error) (Blakcsle e et al.l 12009). This 
cluster is smaller, denser, more dynamically evolved, and 
more regular in shape than the Virgo cluster, and there- 
fore allows us to study the properties of the nuclei of 
galaxies residing in a new and different environment. 

Other papers in the ACSFCS series have described 
the data reduction procedures used in the survey 
(jJordan et al.ll20~07aL hereafter Paper I) , syst ematic vari- 
ation s in the central structure of galaxies (|C6te et al.l 
120071 hereafter Paper II), the logarithmic slope of the 
galax y central surface brightness profiles (Gl ass et al.l 
120111 hereafter Paper IV), and the use o f surface bright- 
ness fluctuations as a distance indicator (Blakcsle e et al.1 
120091 hereafter Paper V). Paper III (2012, in prep.) of 
the ACSFCS will present a detailed isophotal analysis 
of the ACSFCS galaxies, including their dust properties, 
axial ratios, 2D structure, total magnitudes, colors, and 
surface brightness and color profiles. Papers studying the 
properties of globular clusters (GCs) in ACSFCS galaxies 
have examined their half-light radii ([Masters et al.ll2010l 
herea fter Paper VII) , luminosity function (jVillegas et al.l 
120101 hereafte r Paper VIII), color-magnitude relation 
( Mieske et alj 12010), he reafter Paper IX), and color gra- 
dients ( Liu et al.ll2011l hereafter Paper X). 

The outline of this paper is as follows. In f}2] we de- 
scribe the observations and methodologies used to mea- 
sure photometric and structural parameters for the nu- 
clei; in 521 we examine the nucleus properties, including 
their frequency of nucleation, luminosity function, sizes, 
surface brightness parameters, and colors; in f}4] we put 
our results into the context of current formation scenar- 
ios; and in §2] we summarize our main results. An ap- 
pendix presents a comparison of ID and 2D methods for 
measuring photometric and structural parameters of nu- 
clei and their host galaxies. 

2. OBSERVATIONS AND ANALYSIS 

The ACSFCS sample was constructed by selecting 
all galaxies from the FCC with: (1) blue magnitudes 
Bt < 15.5; and (2) early-type morphologies: i.e., E, SO, 
SB0, dE, dE,N or dS0,N. These morphological types were 
taken directly from Ferguson (1989 ) which are in turn 
based on the classification scheme of lSandage fc Binggelil 
(1984). In addition to the 42 FCC galaxies that met these 
criteria, two ellipticals that lie just beyond the FCC sur- 
vey region (NGC 1340 and IC 2006) were added, giving a 
total of 44 targets. Unfortunately, due to a shutter failure 
during execution, no images were obtained for FCC 161 
(NGC 1379). Our final sample therefore consists of 
43 early-type galaxies, which is complete (apart from 
FCC 161) down to a limiting magnitude of Bt ~ 15.5 
mag (Mb ~ —16.0 mag). For all galaxies in this survey, 
membership in the cluster has been confirmed through 
radial velocity measurements. More details on the sam- 
ple can be found in Papers I and III. 

In §4.1[ we will compare our results to a sample of 
galaxies and nuclei from the ACSVCS, which consists 
of 100 early-type members of the Virgo Cluster. That 
survey was magnitude limited down to Bt ~ 12 mag 
(Mb ~ — 19 mag) and 44% complete down to its limiting 
magnitude of Bt ~ 16 mag (Mb ~ —15 mag). Both the 



Fornax and Virgo galaxies were observed with the ACS 
using Wide Field Channel (WFC) mode with the F475W 
and F850LP filters, which correspond closely to the g- 
and z-band filters in the Sloan Digital Sky Survey (SPSS) 
system (see, e g.. iFukugita et alJ ll9M lYork et al.ll2000t: 
ISirianni et al.H2005h . 

Basic data for the ACSFCS galaxies are presented in 
Table [2 The ACSFCS ident ification number, the FCC 
number from IFergusonl (|1989f ) , and any alternate names 
are reported in the first three columns. The table is 
ordered by increasing FCC blue magnitude, Bt, which 
is given in column 4. In calculating absolute magni- 
tudes, we used the individual surface brightness fluctu- 
ation (SBF) distances measured in Paper V. Beginning 
in §3, all reported m agnitudes are extinct ion corrected, 
using dust maps from [Schlegel et al. (1998|), with the ra- 
tios of tota l-to-selective a b sorpt ion in the WFC filters 
taken from ISirianni et al.l (|2005h : the adopted B-band 
extinctions are shown in column 5. The galaxy g- and 
z- band surface brightness at a geometric mean radius 
of 1", measured by spline interpolation, are recorded in 
columns 6 and 7. Note that all HST/ACS magnitudes 
quoted in this paper are AB magnitudes. 

The final two columns in Table [T] g ive the classifica - 
tions of the galaxies as nucleated from IFergusonl |l989), 
and the ones derived from our surface brightness profile 
analysis. The parameterization of these profiles is dis- 
cussed in i j2.ll and the fitting methods used are outlined 
in §2.21 Finally, the nucleus properties obtained from the 
above procedure are described in £12.31 Additional infor- 
mation about our program galaxies, such as coordinates 
and morphological classifications, can be found in Papers 
I and III. 

2.1. Parameterization of the Surface Brightness 
Profiles 

As stellar nuclei, which are the focus of this study, are 
found in the luminous central regions of their host galax- 
ies, accurately modeling the underlying galaxy surface 
brightness is necessary to measure their photometric and 
structural parameters. Indeed, for the faintest nuclei, 
or for some nuclei embedded in high surface brightness 
galaxies with steeply rising brightness profiles, this can 
be important for even identifying a central nuclear com- 
ponent (see Appendix A of[C06). Using t he IRAF task 
ellip se, which is based on the algorithm of lJedrzeiewskil 
(1987), elliptical isophotes with logarithmically increas- 
ing semi-major axis length were fitted to the galaxies. 
In most cases, all ellipse parameters (center, ellipticity, 
and position angle) were allowed to vary. However, to 
achieve convergence, the galaxies with large amounts of 
central dust required the ellipse centers to be held fixed 
throughout the fit (FCC 335, FCC 119, FCC 90), as well 
as the position angles and ellipticities while fitting the 
innermost areas (FCC 119 and FCC 90), where the fixed 
parameter values were determined by ellipse fits to the 
outer regions (R e > 5"). For more detail s on th e fit- 
ting procedures, see §3.2 of iFerrarese et al.1 (|2006bl ) and 
Paper III. 

The results from the ellipse isophotal analysis 
were used to derive azimuthally-averaged radial surface 
brightness profiles, which were then fitted using one of 
three different parameterizations for the global surface 
brightness profile. The first parameterization is the well 
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TABLE 1 
Basic Data for ACSFCS Galaxies 



ID Name Other B T A B Mg(l") /**(!") Class Class 

(mag) (mag) (mag/D") (mag/n") (FCC) (ACS) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 



1 


21 


N1316 


9.06 





090 


15.61 


14.12 


N 


cS 


2 


213 


N1399 


10.04 





056 


16.78 


15.17 


N 


cS 


3 


219 


N1404 


10.96 


() 


049 


16.45 


14.88 


N 


cS 


4 


1340 


E418-G005 


11.23 


(} 


077 


17.00 


15.56 


N 


cS 


5 


167 


N1380 


10.84 


() 


075 


16.88 


15.32 


N 


SI 


6 


276 


N1427 


11.79 


f) 


048 


17.07 


15.59 


N 


SI 


7 


147 


N1374 


11.95 





060 


17.14 


15.58 


N 


SI 


8 


2006 


E359-G007 


12.59 


o 


048 


17.72 


16.18 


N 


S2 


g 


83 


N1351 


12.33 


o 


06 1 


17.35 


15.83 


N 


SI 


10 


184 


N1387 


11.77 


o 


055 


16.70 


15.05 


N 


SI 


11 


63 


N1339 


12.77 





057 


17.13 


15.56 


N 


S2 


12 


193 


N1389 


12.59 





046 


17.34 


15.88 


N 


S2 


13 


170 


N1381 


12.91 


o 


058 


17.12 


15.62 


N 


S2 


14 


153 


11963 


13.55 


o 


062 


18.32 


16.91 


N 


S2 


15 


177 


N1380A 


13.60 


o 


063 


18.83 


17.58 


N 


S2 


16 


47 


N1336 


13.34 





049 


18.50 


17.11 


N 


S2 


17 


43 


11919 


13.82 





062 


19.99 


18.83 


Y 


S2 


IS 


190 


N1380B 


13.79 


o 


074 


19.32 


17.89 


N 


S2 


19 


310 


N1460 


13.68 


o 


047 


19.32 


17.96 


N 


S2 


20 


249 


N1419 


13.61 


o 


056 


17.68 


16.25 


N 


S2 


21 


148 


N1375 


13.39 





063 


18.20 


17.02 


N 


S2 


22 


255 


E358-G50 


13.99 





025 


19.50 


18.26 


Y 


S2 


23 


277 


N1428 


14.01 


o 


044 


18.84 


17.45 


N 


S2 


24 


55 


E358-G06 


14.23 


o 


043 


19.68 


18.41 


Y 


S2 


25 


152 


E358-G25 


14.13 


o 


044 


20.44 


19.25 


N 


SI 


26 


301 


E358-G59 


14.22 


() 


039 


18.61 


17.31 


N 


S2 


27 


335 


E359-G02 


14.90 


o 


063 


20.40 


19.27 


N 


S2 


28 


143 


N1373 


14.19 





001 


18.39 


16.96 


N 


SI 


29 


95 


G87 


15.01 





064 


20.16 


18.83 


N 


S2 


30 


130 


G99 


15.00 





069 


20.73 


19.39 


Y 


S2 


31 


182 


G79 


15.01 





057 


19.61 


18.18 


N 


S2 


32 


204 


E358-G43 


15.33 





045 


20.50 


19.23 


Y 


S2 


33 


119 


G26 


15.44 





060 


21.35 


20.10 


N 


Sl a 


34 


90 


G118 


15.10 





052 


19.55 


18.76 


N 


S2 


35 


26 


E357-G25 


15.26 





067 


19.80 


19.39 


N 


SI 


36 


106 


G47 


15.34 





046 


19.89 


18.62 


Y 


S2 


37 


19 


E301-G08 


15.81 





085 


21.56 


20.49 


Y 


S2 


38 


202 


N1396 


15.50 





057 


20.71 


19.41 


Y 


S2 


39 


324 


E358-G66 


15.83 





042 


22.16 


21.01 


N 


S2 


40 


288 


E358-G56 


15.82 





025 


21.03 


19.85 


Y 


S2 


41 


303 


NG47 


15.74 





046 


21.63 


20.49 


Y 


S2 


42 


203 


E358-G42 


15.82 





051 


21.50 


20.28 


Y 


S2 


43 


100 


G86 


15.75 





062 


22.18 


21.08 


Y 


S2 



Column key: 

(1) ACSFCS Identification number; 

(2) Galaxy name, mainly from the Fornax Cluster Catalog (FCC) of Ferguson (1989); 

(3) Alternative names in the NGC, ESO or IC catalogs; 

(4) Total blu e magnitude from AC SFCS (Paper III); 

(5) A B from lSchlegel et alj fll998D ; 

(6) — (7) g- and z-band surface brightness measured at a geometric radius of 1 ; 

(8) Nuclear classification in the FCC: Y — nucleated, N — non-nucleated; 

(9) Nuclear classification in ACSFCS: cS— core-Sersic (non-nucleated), SI — Sersic (non- 
nucleated), S2— doublc-Scrsic (nucleated) 

a Due to the offset of the nucleus and the amount of central dust, the nucleus parameters 
for FCC 119 were derived using a King profile fit to the ACS image. 



known Sersic profile (jSersid Il968ft , a three parameter 
model which has the form 



Is(R) = i e exp < -b n 



R 



4' 



(I) 



where I e is the intensity at the effective radius, R e , and 
the Sersic index, n, characterizes the overall shape of 
the light profile. The constant b n is defined such that 
r(2n) = 2j(2n, b n ), where T and 7 are the complete an d 
incomplete gamma functions, respectively (|Ciottil fl991). 
For lower values of n, the Sersic profile is shallow in the 
inner regions and steep in the outer regions; n = 1 pro- 
duces a pure exponential profile, which generally pro- 



vides a reasonable fit to dwarf galaxies. Higher values 
of n yield functions which are steep in the inner regions 
and extended at large radii, with a less pronounced radial 
dependence on slope; these profiles generally fit bright el- 
lipticals quite well (i.e., n — 4 reduces to a classical de 
Vaucouleurs profile). 

Historically, these two types of profiles have been used 
to separately parameterize dwarfs and giants. How- 
ever, more complete studies of galaxies have found 
that n actu ally varies steadily with galaxy lumi- 
nosit y fe.g.. [Graham fc Guzmanl 120031: iGavazzi et all 
20051: iFerrarese et all l2006bl; iKormendv et all 120091: 
Misgeld fc Hilkerl 1201 ll McLaughlin et al. 2012, in 
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Fig. 1. — Five galaxies from the ACSFCS chosen to illustrate systematic trends in central and global structure along the luminosity 
function. Top Row: ACS/F475W images showing the inner 15"xl5" (~ 1.5 kpc xl.5 kpc) for each galaxy. Middle Row: Azimuthally- 
avcraged surface brightness profiles and the best-fit PSF-convolved models. For FCC 213, the best-fit "core-Sersic" model is shown. For 
FCC 83, the solid curve shows a fitted Sersic model. The three remaining galaxies — which show central nuclei, or luminosity "excesses" 
relative to an underlying Sersic model — are fitted with double-Sersic models, with the dotted curves showing the separate nuclear and 
global components. The arrow in each panel is drawn at 2% of the effective radius of the galaxy {Paper IT) . Bottom Row: Deprojected 
lumi nosity profi les (which represent the true 3-dimcnsional density distribution without any PSF convolution) for the same five galaxies 
from lPaper IVI . The solid curves show the deprojected profiles corresponding to the solid curves shown in the middle row. Dashed curves 
show the profiles corresponding to the inward extrapolation of the Sersic models that best fit the outer (galaxy) profile. 



prep.). In what follows, we will refer to these single- 
component parameterizations as SI models. 

Although the Sersic profile describes the outer compo- 
nent (typically beyond a few percent of the the effective 
radius) of galaxies remarkably well — a consequence of 
the wide range in concentration, spatial scale and sur- 
face brightness that is possible by varying n, R e and 
I e , respectively — there can be variations in the cen- 
tral structure that cannot be accounte d for in th is simple 
model (see, e.g., Figures 1 and 2 of lPaper"TH ) . Specifi- 
cally, the brightest ellipticals tend to show a luminos- 
ity deficit in their central regions; for these objects, the 
six-parameter "core-Sersic" model ()Graham et al.ll20"03() 
provides a good description of their surface brightness 
profiles. The core-Sersic model, referred to hereafter as 
a cS profile, can be written as 



Ics(R) = I' 



where 





■y/a 




j^a \ I /an 




exp 


[-H R 


? ) 








(2) 


I h 2~ lla exp 


bn (2 1 / 


\ 1/n ' 
Rb/Re ) 


■ 0) 



This parameterization consists of the usual Sersic profile, 
with effective radius R e and Sersic index n, outside of a 
"break" radius Rb (where the intensity is /&). At Ri,, the 
outer profile transitions to an inner power-law component 
with slope 7, according to the "sharpness" parameter a 
(where smaller values translate to smoother transitions) . 

By contrast, most of the low- and intermediate- 
luminosity galaxies in our sample show evidence for a 
luminosity excess in their cores which is, by definition, 
the signature of a central nucleus (see Appendix A of 
ICote et all 12006) Fl A central excess in the surface 
brightness profile can then be modeled by adding a sec- 
ond Sersic component. This double-Sersic profile (which 
we denote hereafter as an S2 profile) has the form 



Is2(R) = h,i exp < -b„,i 
+ Z e ,2 exp <^ -b n o 



(A) 



l/ni 



l/«2 



(4) 



13 The ACSVCS finding of "luminosity excesses" in Virgo clus- 
ter galaxies relative to the inward extrapolation of Sersic mod- 
els fitted to the outer profiles was subsequently confirmed by 
IKormendv et aEI (2009) who reanalyzed a subset of the ACSVCS 
sample. 
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Fig. 2. — F475W (g-band) images of the inner 10" X 10" (~ lkpc X lkpc) regions of the ACSFCS galaxies. The galaxies are arranged 
in order of increasing blue magnitude (i.e., decreasing luminosity) from left to right, and from top to bottom. Each galaxy's FCC number 
is displayed in the top left, and the bottom left denotes the model used to fit the galaxy, either SI (Sersic), cS (core-Sersic), or S2 
(double-Sersic). 



where the enumerated subscripts indicate the Sersic pa- 
rameters for the outer and in ner c omponents. 

It should be noted that, in lC06l double-Sersic profiles 
were not used to fit the nucleated galaxies. Inste ad, the 
central nuclei we re represented by King profiles (|Michiei 
119631: iKmsl I1966D . while the outer component was rep- 
resented by either a core-Sersic or Sersic profile. Our 
decision to use a double-Sersic parameterization in the 
ACSFCS analysis is motivated by two considerations. 
Firstly, modeling the inner component with the Sersic 
profiles allows for a diversity of possible physical systems, 
due to the range of the Sersic parameter (see above). 
For n ~ 1, the profile is a pure exponential and is thus 
suitable for embedded disks, whereas n ~ 2 represents 
Galactic GCs quite accurately, and presumably, nuclear 
star clusters as wel l . This is supported by the findings of 
iGraham fc Spitleil (|2009f) . who measured Sersic indices 
of of n = 3.0, 2.3, and 1.6 for the nuclear star clusters of 



the Milky Way, M 32, and NGC 205, respectively; and by 
ISeth et all (|2010f) who observed a Sersic index of n ~ 2 
for NGC 404. Secondly, the use of Sersic profile for both 
the inner and outer components allows straightforward 
and convenient comparisons of their respective structural 
properties. 

The overall trends described here are illustrated in Fig- 
ure [TJ which shows the systematic variations in the core 
and global structure of ea rly-type galaxies along the l u- 
minosity function ( see also ICOalFerrarese et al.|[2006bllal ; 



iPaper iHlPaper IVf ). The upper panels in this figure show 
15" x 15" images centered on five representative galax- 
ies from the ACSFCS, arranged in order of decreasing 
luminosity. The middle rows show model fits to the g- 
band surface brightness profiles as described above: i.e., 
FCC 213 (cS), FCC 83 (SI), FCC 277 (S2), FCC 136 (S2) 
and FCC 303 (S2). Note the systematic decline in galaxy 
surface brightness from left to right, and the emergence 
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of an increasingly prominent central nuclear component 
as galaxy luminosity decreases. At low and intermedi- 
ate luminosities, these luminosity "excesses" (i.e., nuclei) 
relative to the underlying galaxy model correspond to 
a steady steepening of the three-dimensional luminosity 
density on sm all scales, as shown in the lower panel of 
Figure Hlffrom lPaper IVT ) . Images and brightness profiles 
for the full sample of ACSFCS galaxies will be discussed 
below. 

2.2. Fitting Procedure 

As described in iPaper 1 the ACSFCS uses the Lanc- 
zosS kernel for drizzling rather than the Gaussian kernel 
which was selected for the ACSVCS. Due to the slightly 
larger distance of the Fornax cl uster — 20.0 vs. 16.5 
Mpc ijMei et all 120071: IPaper VP — and the fact that 
some of the V irgo nuclei were only margin ally resolved in 
the ACSVCS ([CM iFerrarese et al.ll2006bllal ). the sharper 
point spread function (PSF) possible with the LanczosS 
kernel was deemed to be more important for the ACS- 



FCS galaxies than the Gaussian kernel's ability to repair 
bad pixels. 

New PSFs for the ACSFCS were constructed in an 
identical manner using more than a thousand stars from 
the GO-10048 and GO-10375 programs to obtain photo- 
metric calibrations of the Galactic GC 47 Tucanae (PI = 
J. Mack). Using multiple observations allowed PSFs to 
be extracted from data that were acquired no more than 
two months away from the ACSFCS observation times; 
this proved to be important since on 2004 December 20, 
the HST secondary mirror wa s moved by 4.6 imq. . 

After running KINGPHOT ([Jordan et alj|2005[) on the 
GC candidates identified in the ACSFCS images F^l it was 
found that, for a subset of galaxies (FCC 213, IC 2006, 
FCC 193, FCC 249, FCC 277, FCC 19, and FCC 202), 
the mean half-light radius for GC candidates was sig- 
nificantly larger in the g-band than in the z-band: i.e., 



14 KINGPHOT fits two-dimensional, PSF-convolved King mod- 
els to candidate GCs in the ACS images. 
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Fig. 2. — Continued 

by roughly 0.5 pixels in F475W, which is much larger 
than the < 0.1 pixel diff erences found in the ACSVCS. 
lAnderson feKmel (|20"06l ) showed that the WFC PSF ex- 
hibits unpredictable variations on orbital timescales, par- 
ticularly in the bluer filters, with differences in flux val- 
ues of up to ~ 10% in the central regions. To correct the 
seven galaxies whose imaging suffered from this variabil- 
ity, stellar sources in the individual images were used to 
adjust empirically the 47 Tucanae P SFs. Full details on 
this procedure are given in lPaper 1 

The azimuthally-averaged, one-dimensional (ID) sur- 
face brightness profiles were fitted using a \ 2 minimiza- 
tion scheme to determine if a Sersic or core-Sersic model 
was most appropriate. If visual inspection of the im- 
ages and/or surface brightness profiles revealed a nucleus, 
then an S2 parameterization was adopted. At each it- 
eration of the fitting procedure, the models used were 
convolved with the PSF in two dimensions (assuming 
spherical symmetry) , and both the models and PSF were 
oversampled by a factor of 10 with respect to the ACS 
pixel size (i.e. they were sampled every 0'.'005). 

All profile parameters, except for intensity, were first 
fitted to both bandpasses simultaneously. These prelimi- 
nary values were then used as initial guesses for the inde- 
pendent g- and z-band fits for most of the galaxies, with 
the exception of those with high central surface bright- 
ness that appear to be nucleated. In these galaxies, the 
nuclei are often quite extended, and difficult to differen- 
tiate from the underlying galaxy light; thus, only the in- 
tensity parameters were allowed to vary between the two 



bands. As many previous investigators have noted, it is 
possible to measure reliably the total magnitudes and ef- 
fective radii of marginally resolved stellar systems (i.e., 
star clusters, nuclei) using HST imaging, whereas the 
concentrations can us ually be constrained with consider- 
ably l ower precision (iKundu fc Whitmorelll998t iLarsenl 
1199a iCarlson fc Holtzmanl 120011; Uordan et all 120051 7 
This is understandable given that the measurement of 
concentration (or Sersic index) for a stellar system re- 
quires the curvature of the profile to be measured on 
scales smaller than the PSF. Fortunately, the derived 
radii and magnitudes are quite insensitive to Sersic in- 
dex, at least insofar as the adopted model is an accurate 
representation of the actual nuclear profile. 

A co nservative resolution limit of 0'.'025 was estimated 
in lC06l based on the half-light radii of King models fit to 
stars classified as unresolved by KINGPHOT, and from 
the size of the cen tral n on-thermal point source found in 
VCC 1316 (M87). [C0l further showed that most of their 
detected nuclei were more extended than point sources, 
by fitting point source profiles in addition to King pro- 
files, and comparing the residuals. Four of the nucleated 
galaxies (FCC 301, FCC 249, FCC 255, FCC 95) in Ta- 
ble [2] have best-fit effective radii that are measured to be 
smaller than our resolution limit in one, or both, pho- 
tometric bands; these nuclei are thus unresolved — or 
nearly so — in our HST imaging. 

After some experimentation, we have estimated uncer- 
tainties on the fitted parameters for the nuclei (and their 
host galaxies) using a Monte Carlo approach in which 
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Fig. 3. — Azimuthally-averaged surface brightness profiles for the ACSFCS galaxies. As in Figure [2] galaxies are ordered by blue 
magnitude, with luminosity decreasing from left to right, and from top to bottom. The black points plot the measured g-band profiles. 
The red curves show the fitted models with the two separate components (nucleus and galaxy) indicated by the dashed curves; their sum 
is shown as the solid curve. The dotted vertical lines are drawn at a radius of 0.02R e in all cases. The top right label denotes the galaxy 
FCC number, and the three types of fitted models are denoted in the bottom left, by SI (Sersic), cS (core-Sersic), or S2 (double-Sersic). 
Note that FCC 167 contains a prominent central dust disk (Figure [2] Paper III), so the models were fitted outside R = 5". 



the g- and z-band surface brightness profiles for each 
galaxy are independently simulated 200 times. We in- 
cluded an amount of noise at each data point in the pro- 
file assuming a Gaussian distribution of errors and using 
the uncertainty on the intensity at each point computed 
by ELLIPSE. An additional source of error for the profiles 
comes from the determination of the background level, 
which we have also included by assuming a 10% error 
in the adopted background for each galaxy (estimated 
roughly by the galaxy-to-galaxy sca tter in th e measured 
background levels; see Figure 5 of iPaper 1 ). Errors on 
the magnitudes, colors and effective radii estimated from 
these Monte Carlo simulations are given in Table 1. We 
hasten to point out that these errors do not include pos- 
sible sources of systematic errors, such as errors in the 
PSF, and that they are therefore best viewed in a relative 



sense, and as lower limits on the true errors. 

Finally, as a check on the (ID) method, we also fit- 
ted surface brightness profiles to our galaxies using 2- 
dimensional (2D) techniques. The full results of this 
analysis are described in In brief, the structural pa- 
rameters obtained from the two procedures are largely in 
agreement, but due to the increased difficulty of charac- 
terizing complex structures using 2D fitting, we proceed 
with results from the ID method, which we consider most 
appropriate for this study. 

2.3. Identification of the Nuclei 

The classification of a galaxy as nucleated or non- 
nucleated was performed in the following way. The pro- 
gram galaxies were all fitted with pure Sersic profiles out- 
side of a geometric mean radius of ft' 5 (~ 50 pc). The 
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geometric mean radius was derived from the fitted ellip- 
tical isophotes, and is thus defined as R = a(l — e) 1 ^ 2 , 
where a is the semi-major axis, and e is the ellipticity. If 
an inward extrapolation of this profile revealed an excess 
of light in the center, then the full profile was refitted by 
adding a second Sersic component, and the galaxy was 
thus considered nucleated and classified as S2. In gen- 
eral, the level of nucleation was slightly greater in the 
g-band, as the nuclei are often found to be somewhat 
bluer than their hosts (see §3.6|) . 

One of our program galaxies, FCC 119, appears to have 
a distinct nucleus offset that is from its photoccntcr by 
~ 0'.'7. Due to the presence of dust in the inner re- 
gions of the galaxy, the ellipse centers were held fixed 
to the photocenter throughout the fit; thus, the nucleus 
is not apparent in the ID surface brightness profile (dis- 
cussed below) . We therefore use parameters derived from 
a KINGPHOT fit to this object, and consider this galaxy 
to be nucleated for the remainder of our analysis. 

Galaxy classifications as nucleated or non-nucleated in 



the FCC, and our revised classification, are presented in 
columns 8 and 9 of Table [T] In Table [3J we record the 
parameters of the Sersic profile fit to the nucleus of all 
S2 galaxies, as well as the KINGPHOT fit to FCC 119. 
Specifically, we have measured the g- and z-band inte- 
grated nucleus magnitudes (columns 3 and 4), integrated 
and 4-pixel radius aperture nucleus colors (columns 5 and 
6), and g- and z-band nucleus half-light radii (columns 7 
and 8). Error estimates for each of these parameters are 
also included in this table, derived using the Monte Carlo 
approach described in £12.21 We have also calculated, by 
integrating the Sersic profiles, the fraction of luminosity 
occurring outward of B > 0'.'5 (columns 7 and 8). 

Although the nuclei half-light radii in a few galaxies 
were measured to be somewhat larger in the g-band than 
in the z-band (i.e. FCC 310, FCC 177, FCC 95), we 
note that these are not the same galaxies that suffered 
from the variable PSF discussed in $2.21 for the most 
part, these differences reflect the fact that size measure- 
ments are particularly challenging for under-luminous or 
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extended nuclei in galaxies with steeply rising surface 
brightness profiles. For the reasons discussed in £)2.21 
we do not report the best-fit Sersic indices in Table [2 
although we note that the indices for all nuclei in our 
sample have 0.5 < n < 4, with a median of n = 2.0 ± 0.7. 

F475W images for the central 10" x 10" region of the 
program galaxies, where a distinct nuclear component 
is often discernible, are shown in Figure [5J The FCC 
number of the galaxy is labeled in each of the panels, 
along with the type of profile fitted; S2 therefore indi- 
cates that the galaxy was considered to be nucleated 
(that is, fitted with a double-Sersic profile). Individual 
fits to the azimuthally-averaged (/-band surface bright- 
ness profiles are shown in Figure [3] These im ages illus- 
trate the systematic trend noted in iPaper III , in which 
the central regions of early- type galaxies transition from 
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to a two-component structure (nucleus+galaxy) as one 
moves down the luminosity function: i.e., toward fainter, 
and lower surface brightness, galaxies. 

3. RESULTS 

In the following section, we analyze the properties of 
the ACSFCS nuclei derived from the above parameteri- 
zation. Specifically, we examine the frequency of nucle- 
ation (§ 13. ip . offset of the nuclei from their hosts (§ 13. 2|) . 
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Fig. 4. — Top: Luminosity distribution of the 43 ACSFCS pro- 
gram galaxies (open histogram). The overlaid hatched histogram 
shows the distribution of the 31 galaxies classified as nucleated in 
this study, while the solid histogram shows the distribution of the 
12 nucleated galaxies according to the FCC. Bottom: The per- 
centage of nucleated galaxies (/ n ) in this study (open squares) and 
in the FCC (solid squares). The approximate luminosity regimes 
where galaxies show central surface b rightness "deficits" a nd "ex- 
cesses" (i.e., nuclei) are indicated (see Paper II; Pape rTVf) . 



nucleus-to-galaxy luminosity ratio (§ 13 . 3[) . and nucleus 
luminosity function (§ 13 .4|) . Additionally, we compare 
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Fig. 5. — Galaxy surface brightness in the g-band (top) and z- 
band (bottom) measured at a mean radius of 1", plotted against the 
magnitude of the nucleus. The filled circles show the 31 galaxies 
found to be nucleated by this study, while the open squares show 
the 12 galaxies classified as nucleated in the FCC. 



their structural properties and scaling relations f? 13. 5p , 
as well as colors (^ I3.6[) , with those of their host galaxies 
and GCs. 

3.1. Frequency of Nucleation 

Only 12 out of our 43 program galaxies were classified 
as nucleated in the FCC, which sets the frequency of 



nucleation at /„ w 28%03 Column 8 of Table Q] shows 
the classification as nucleated or non-nucleated in the 
FCC. These can be compared to our classification in the 
ACSFCS, where the use of the double-Sersic (S2) model 
indicates that we consider the galaxy to be nucleated. 
We find all galaxies previously classified as nucleated in 
the FCC to be nucleated in our sample, as well as an 
additional 19 objects, for a total of 31/43 galaxies, or 
/„ « 72%. 

The cause of this sharp rise in frequency of nucleation 
can be attributed to both observational selection effects 
and the definition in our analysis of a nucleus as a central 
luminosity excess relative to a fitted galaxy model. In the 
top panel of Figure [H the open histogram shows the lu- 
minosity of all of the program galaxies, while the hatched 
and solid histograms denote those found to be nucleated 
in the ACSFCS and the FCC, respectively. The bottom 
panel of Figure 2] plots /„ as a function of luminosity 
for the two surveys. The ACSFCS uncovers many more 
nuclei in more luminous host galaxies, as the high resolu- 
tion of the WFC allows us to resolve nuclei in their high 
surface brightness cores. This selection effect is explored 
further in Figure O 

Galaxy surface brightness at a geometric mean radius 
of R = 1" (sa 97 pc) was calculated using linear spline 
interpolation, in the g- and z-bands. By measuring sur- 
face brightness at a constant radius (rather than at some 
function of the effective radius), the result is a model- 
independent measure of central surface brightness, at a 



15 We include NGC 1340 and IC 2006 in this calculation; al- 
though they do not appear in the catalog of Ferguson (1989), both 
have "E" classifications in NED (i.e., non- nucleated ellipticals). 
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Fig. 6. — Top: Projected offset between the nucleus and the 
galaxy photocenter in the g-band, plotted against host galaxy mag- 
nitude. Offsets were calculated using our ID (black circles) and 2D 
(blue triangles) analyses. The two dotted red lines show offsets one 
and ten ACS/WFC pixels (0'.'05 and Of.' 5). The black short-dashed 
line and blue long-dashed line represent the best-fit relation for ID 
and 2D offsets, respectively. Bottom: Same as above, but for the 
z-band. 

distance large enough to avoid the contribution from a 
typical nucleus, if present. Figure [5] plots the integrated 
nucleus magnitude derived from the S2 fit against galaxy 
surface brightness measured at a distance of 1". The 
filled circles and open squares show the measurements 
for galaxies classified as nucleated in the ACSFCS and 
FCC, respectively. Clearly, the nuclei that went unde- 
tected in the earlier (photographic) survey come in two 
forms: bright nuclei that are embedded at the centers of 
galaxies with intermediate luminosity (which also have 
steeply rising profiles, see Figure [3]) and faint nuclei be- 
longing to the lowest luminosity galaxies. Needless to 
say, it is possible that we too may be missing some nu- 
clei, so we take /„ « 72% as a lower limit on the true 
frequency of nucleation in the ACSFCS sample. 

3.2. Offset Nuclei 

The offset of each nucleus from its host galaxy photo- 
center was measured for 28 of our 31 nucleated galax- 
ies. For FCC 335, FCC 119, and FCC 90, the ellipti- 
cal isophote fitting was performed with the ellipse cen- 
ters held fixed, as convergence could not be otherwise 
achieved due to dust in their central regions. Thus, off- 
sets for these nuclei could not be measured using the 
technique described below, although we did examine 
their offsets using our 2D (GALFIT) analysis, and they 
are also included in this section. 

For the remaining 28 galaxies in qu estion, an analysis 
sligh tly different to that used in IC06I was performed. In 
IC061 the galaxy photocenter was determined by taking 
the mean of the positions of all fitted isophotes satisfy- 
ing 1" < R < R e . However, because of the possibility 
that isophotes might drift from the center due to bright 
sources in the field of view (causing an artificial offset in 
the photocenter calculated using the above method), we 
have adopted a different procedure in this work, where 



the photocenter and its error was determined by running 
ellipse to fit a single isophote with a sem i- major axis 
length of approximately R e /2. As in IC061 the position 
and error of the centroid of the nucleus were taken as the 
smallest fitted ellipse from the full ellipse run. We note 
that the geometrical parameter errors output by ellipse 
are calculated from the errors of the harmonic fit, with 
the first and second harmonics removed. The errors from 
the photocenter and centroid were then added in quadra- 
ture to obtain the total error on the offset. The results 
of this procedure are plotted as the black filled circles in 
Figure [6] We find that almost all of the galaxies in our 
sample have an offset of less than Of! I. The four that do 
have a larger offset (FCC 63, FCC 193, FCC 177, and 
FCC 277) have AR n > O'.'l in only one of the two bands. 

We generally observe the offsets from our 2D analysis 
to be larger than those determined using our ID method. 
This is due to the fact that the 2D fitting procedure does 
not allow the ellipse parameters to vary with radius, and 
returns the model that best fits the average parameters 
of the entire galaxy, giving more weight to the outer re- 
gions in the determination of the photocenter. Thus, 
for our 2D analysis, we are not concerned with offsets 
larger than ~ 0'.'5, and only three galaxies are found to 
have offsets larger than this — FCC 119 (0"65), FCC 324 
(0'.'70), and FCC 288 (0"62). FCC 119 is fairly irregular 
in structure, with a significant amount of dust in its core. 
FCC 324 and FCC 288 are both low surface brightness, 
highly flattened galaxies, with no obvious clusters near 
the photocenter that may have caused source confusion 
with what we consider to be the nucleus. We conclude 
that in our sample, at most 10% of the nuclei are offset 
at the level of 0'.'5 or more, consistent with the findings 
of IC06I for the Virgo cluster. 

To measure any trend between offset and galaxy lumi- 
nosity, we perform a weighted least-squares fit to the data 
from Figure [Sj Using the offsets from our ID analysis, 
we obtain 

logAR g = (-0.057 ± 0.070) B T - (0.87 ± 1.06) 
logAi? 2 = (0.19 ± 0.09) Bt — (4.9 ± 1.3), [b) 

and from our 2D analysis, 

logAi? 2D = (0.21 ± 0.18) B T - (4.1 ±2.7). (6) 

The slopes of these relations do not indicate any signif- 
icant trend between offset and galaxy luminosity. The 
errors on the fitted parameters are the standard errors. 

Finally, some of the galaxies in our sample that we do 
not find to be nucleated may, in fact, be "dlrr/dE tran- 
sition" objects, where a nucleus could be in the process 
of formation^ In particular, FCC 152 and FCC 26 are 
irregular in shape and contain many star clusters and 
significant amounts of dust in their central regions. It 
is possible that one or more of these clusters could be 
nucleus progenitors that will migrate inwards through 
dynamical friction (see £)4.3.1|) . 

3.3. The Nucleus-to- Galaxy Luminosity Ratio 

Previous studies of early- type dwarfs (iLotz et al.ll2004J : 
IGrant et al.l 120051 : 1 Graham fc Guzmanl 120031 ). including 

16 A prototype for this class is VCC1512 in the Virgo cluster 
which contains a prominent central excess that is composed of blue, 
densely packed star clusters. 
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Fig. 7. — Top: Nucleus magnitude plotted against host galaxy 
magnitude for nucleated galaxies in the ACSFCS; results for the 
g and z bands are shown on the left and right, respectively. The 
lines show the best-fit relations, with the slope held fixed at unity 
(solid) and allowed to vary (dashed). Bottom: Nucleus-to-galaxy 
luminosity ratio, r], against host galaxy magnitude for the g-band 
(left) and z-band (right). The solid and dotted lines show the mean 
value of r\ and its ±lcr limit, while the dashed line shows the best 
fit relation given by the dashed line in the upper panel, recast in 
terms of log(r;) and host magnitude. 



TABLE 3 

Nucleus-to-Galaxy Luminosity Ratios 

Sample Band a 1 ^ 1 2 < lo g V) a 
(mag) (mag) (dcx) (dex) 

ACSFCS g 0.90 ± 0.17 7.27 ± 2.45 5.79 ±0.15 -2.31 0.32 

ACSFCS z 1.07 ± 0.16 5.37 ± 2.12 6.20 ± 0.16 -2.49 0.35 



IC061 found that nucleus brightness increases with host 
galaxy brightness. Similar relations are known to exist 
for the nuclear clusters in late-type galaxies (see, e.g., 
ICarollo et aL1[l998l: iBoker et alj|200l . A plot of nucleus 
versus host galaxy magnitude, the latter calculated by 
integrating the Sersic profile best-fitting the main galaxy 
component over all radii, is shown at the top of Figure [JJ 
Weighted best-fit linear relations of the form: 

3nuc = aggai + P 

z nuc = azg a i + fi (7) 

were fitted to the data, where g n and z n are nuclei mag- 
nitudes and <7 g and z g are the galaxy magnitudes. The 
best-fit parameters (ai, fa, fa) are given in Table 02 
where the quoted errors are the standard errors. Results 
are given for the two cases of fixing the slope at = 1, 
and allowing it to vary freely (shown respectively as the 
solid and dashed lines in the upper panel of Figure [7]). 

Since the best-fit slope of the nucleus-galaxy luminos- 
ity relation is very nearly one, we consider the possi- 
bility of a constant nucleus-to-galaxy luminosity ratio, 
T) = C n /C g , where C n and C g are nucleus and galaxy 
luminosity, respectively. In the bottom of Figure [7j 77 
is plotted as a function of host galaxy magnitude in the 
same band. The values of the weighted means and stan- 
dard deviations are given in Table |3l while the weighted 
mean ratio and standard error on the mean from both 



bands is 

(77) = 0.41% ± 0.04%. (8) 

This is 0.11 ± 0.06% larger than the value of (77) = 
0.30 ±0.04% found m[cM(a 1.9a discrepancy). At first 
glance, this might suggest that, at a given luminosity, 
early-type galaxies in Fornax were slightly more efficient 
in assembling their nuclei than those in Virgo; however, 
the difference is due to the use of Sersic rather than King 
models in fitting the ACSFCS nuclei. As Sersic profiles 
with even moderate n have somewhat extended wings, 
they increase the inferred luminosity of the nuclei rela- 
tive to the King models (whose defining characteristic is 
a tidal truncation radius) . Re- fitting the Virgo data with 
S2 profiles confirms this conclusion — in § 14.1.21 where 
the new fits to Virgo are presented, we find agreement 
between rj for both clusters. 

Finally, due to the definition of ij, the best-fit relation 
from Equation [7J can be recast in terms of log(?7) and 
galaxy magnitude, where a v = —0.4 (a\ — 1) and fa — 
—OAfa. This relation is plotted as the dashed line in 
the bottom panels of Figure [Jj Although we do not see 
any trend between 77 and galaxy luminosity, we note that 
in a study o f galaxies containing bot h a nucleus and a 
black hole bv lGraham fc Spitlerl (|2009f) , it was found that 
the ratio of total CMO mass to spheroid mass tended to 
decrease in more massive galaxies. 

3.4. Luminosity Function 

One mechanism for the formation of galaxy nuclei 
is through multiple mergers of GCs t hat sink to the 
galaxy center by d y namic al friction (e.g..lTremaineill976l: 
Capuzzo-Dolcettal 119931: iCapuzzo-Dolcetta fc Tesseril 
1999HLotz et al.ll2001b iBekki et al.ll2004t l. A comparison 
of the luminosity function of our nuclei with that of the 
GCs identified in the ACSFCS can offer some insight 
into this process. In Figure [SI we present the results 
of a weighted maximum-likelihood fit to the luminosity 
functions of the nuclei, using a normalized Gaussian: 

$(m° uc ) cx exp [-« c - fh° nuc )/2al uc ] . (9) 

While this choice of parameterization is commonly used 
for GCs, there is no physical reason that the nuclei should 
have a Gaussian distribution. It is, nevertheless, a use- 
ful departure point for the purpose of comparison with 
the GCs. To parameterize the GC luminosity function, 
we also performed a maximum-likelihood fit of a normal- 
ized Gaussian, using the GC turnover magnitudes for 
each galaxy (whic h have bee n corrected for complete- 
ness), taken from lPaper Villi . Each turnover magnitude 
was weighted by the number of GCs in the galaxy. 

Our GC sample consists of ~ 2000 candid ates with 
probability index V gc > 0.5 (see iPaper Villi for more 
details on the GC probability index, and a detailed study 
of the GC luminosity functions). The best-fit parameters 
for both nuclei and GCs are given in Table [H where the 
errors on the fitted parameters are the standard errors. 
We find the luminosity function of the nuclei to be both 
brighter, and have a greater spread, than that of the GCs. 
The difference in the means is A = 4.03 mag in each 
band; that is, on average the nuclei are ~ 40 x brighter 
than a typical GC. 

In reality, since we find the nucleus-to-galaxy luminos- 
ity ratio to be roughly constant, the nucleus luminosity 
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Fig. 8. — Luminosity function for the nuclei (closed black circles) 
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plotted for comparison. Both luminosity functions are derived by 
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magnitude was weighted by the number of GCs in the galaxy. 




TABLE 4 

Nucleus and Globular Cluster Luminosity Functions 



Sample 


Band 
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Fig. 9. — Top left: A Schechter function plotted over the B- 
band galaxy luminosity distribution for all early-type galaxies in 
the FCC (open circles) and the ACSFCS sample (filled circles). Top 
right: The previous Schechter function truncated at Bt = 12.5 and 
16, and reduced by 91%, so that it represents the nucleated galax- 
ies in our sample. Bottom left: The previous Schechter function 
shifted by +6.0 mag, corresponding to a constant (rf) of 0.41%, and 
—0.4 mag to convert from B to g. It should now be roughly cor- 
respond to the nucleus luminosity distribution, although without 
taking into account the scatter. Bottom right: The luminosity dis- 
tribution of the nuclei (solid line), the previous Schechter function 
(dashed line), and the same Schechter function convolved with a 
Gaussian of <r = 0.87 mag, the unweighted standard deviation of 
rj g (dotted line). 



Fig. 10. — Distribution of half-light radii of the 31 nuclei identi- 
fied in this study (black filled circles), as well as the candidate ACS- 
FCS GCs (blue open triangles), measured in the the g- (top) and 
z-bands (bottom). The vertical dotted lines indicate the adopted 
resolution limit of ~ 0'.'025. 

function should reflect that of the host galaxies (albeit 
with more scatter), and is most likely parameterized by 
a Schechter function truncated on both ends — on the 
bright end because we find no bright galaxies that are nu- 
cleated, and on the faint end because our sample is mag- 
nitude limited. To illustrate this, we show in Figure [9] a 
Schechter function overlaid on the .B-band galaxy lumi- 
nosity distribution of the FCC and the ACSFCS. We then 
apply cutoffs at Bt = 12.5 and 16 mag (so that we are 
left with the magnitude range of the nucleated galaxies in 
our sample), and scale it down by 91% (since we have 34 
galaxies in this magnitude range, 31 of which are nucle- 
ated), and shift the Schechter function over by +6.0 mag 
(which corresponds to an (rj) of 0.41%) and —0.4 mag 
(to convert roughly from B to <?). Finally, we convolve 
it with a Gaussian with a = 0.87 mag, the unweighted 
standard deviation of 77, to take the scatter around {if) 
into account Q The resulting function plotted over the 
nucleus luminosities shows good agreement, apart from a 
few bright outliers; specifically, the nuclei from FCC 63, 
FCC 47, and FCC 148 have B T < 17. Each of these 
galaxies have a complex central structure that may be 
contributing to their brightness, either by causing accu- 
rate nucleus parameterization to be more difficult, or due 
to the fact that, in these cases, the nucleus themselves 
may be more complex. 

3.5. Structural Properties and Scaling Relations 

17 We use the unweighted standard deviation in this case, be- 
cause we are looking to reproduce the observed rather than intrinsic 
scatter. 
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Fig. 11. — Scaling relations of nuclei compared to galaxies and 
globular clusters from the ACSFCS and ACSVCS surveys (with 
data for galaxies from McLaughlin et al. 2012, in prep.). Top: 
Stellar mass surface density plotted against stellar mass. The 
dotted line shows the relation for GCs, which have R e ~ 3 pc 
llJordan ct al. 2005; Paper VIII). The dashed line shows the re- 
lation calculated from the IBekki et aT] <200ll) finding that nuclei 
assembled from repeated GC mergers have R e oc L o ss . Bottom: 
Effective radius plotted against stellar mass for the same stellar sys- 
tems. The curves are the same as in the previous panel. The solid 
lines show scaling relations of the form R e oc M 5 and R e oc M. 
See text for details. 



Nuclei at the distance of the Fornax cluster are al- 
most never resolved in ground-based imaging, as 1" cor- 
responds to «100 pc at a distance of 20 Mpc. However, 
with ACS resolution it is possible to m easure sizes for 
nuclei as small as R e ~ Of! 025 (see EM). In Figure [H 
we present a comparison of the effective radii of the ACS- 
FCS nuclei and GC candidates. On average, the nuclei 
are larger in size, and have a much greater spread, than 
the GCs, although the considerable overlap between the 
two distributions shows that the most com pact nuclei are 
very nearly the same size as typical GCs ((Paner Vim . 

The two most prominent outliers are the nuclei of 
FCC 63 and FCC 47; both of these nucl ei we re also found 
to be the brightest in our sample (see § !3.4p . Regardless, 
Figure [10] clearly demonstrates that the nuclei have a 
size distribution that peaks at compact sizes, and an ex- 
tended tail populated by larger nuclei. The median sizes 
of the full sample are found to be 0'/073 (7.2 pc) and 
0'.'071 (7.0 pc) in the g- and z-bands, respectively. 

Figure[TT]shows scaling relations for the nuclei from the 
ACSFCS and ACSVCS. In the upper panel, we plot the 
effective mass surface density, S e = A^*/27ri?g, against 
total stellar mass, A^*, calcula ted from the ob served 
(g — z) colors and the relations of I Bell et al.l (|2003). The 
lower panel shows effective radius as a function of stel- 
lar mass. In both panels, we also plot ACSFCS GCs, 
and the sample of early- type galaxies in Virgo, Fornax 
and the Local Grou p from McLaughlin et al. (2012, in 
prep.). As found bv Uordan et al.l (|2005l ). the GCs have 
a size of R e ~ 3 pc that is nearly independent of mass, 
while the early-type galaxies show a smoothly varying 
A4*-R e relation, a reflection of the fact that galaxies 
form a non-homologous family (McLaughlin et al. 2012, 



in prep.). 

This figure highlights several other interesting proper- 
ties of the nuclei. First, there is an obvious similarity 
in the scaling relations of the Fornax and Virgo nuclei; 
we shall return to this point and its implications for nu- 
cleus formation models in §4.11 and 8 34.31 The addition 
of th e ACSFCS nuclei also reaffirms the trend noted by 
IC06I that the nuclei, unlike GCs, obey a size-mass re- 
lation that merges with the GC sequence at low mass. 
For reference, the dashed line in the lower panel of Fig- 
urefTTI shows the predicted scaling relation for nuclei that 
are assembled from repeated GCs mergers, R e oc A4®- 3S 
(|Bekki et al.ll2004D . The corresponding E e -M* relation 
is shown in the upper panel. Based on structural param- 
eters alone, we conclude that the GC merger model is 
broadly consistent with the data (although the extremely 
red colors of the brightest nuclei pose a challenge to this 
model in its simplest form). The two solid lines in the 
lower panel show relations of the form R e oc M.® 5 and 
R e oc M. * , which will be discussed in H4.3.1I 

3.6. Nuclei Colors 

As in [COl, we find a relationship between nuclei col- 
ors and magnitudes with the brighter nuclei having red- 
der colors and residing in more luminous hosts. This is 
shown in Figure 1121 where we plot the nucleus 4-pixel 
aperture colors against g-band magnitudes. Significant 
scatter is seen for the brighter galaxies, which are la- 
beled with their FCC number. This scatter was also 
seen in IC06I . although in Virgo galaxies bright nuclei ap- 
peared to be preferentially red, while in the case of For- 
nax, bright nuclei are seen to scatter to both red and blue 
colors. The increased scatter in the color of the bright nu- 
clei may simply reflect the more complex formation and 
enrichment histories in their inner regions of brighter, 
more massive galaxies: mergers, gas inflow, star forma- 
tion, and GC accretion would naturally lead to a greater 
degree of scatter in the general color-magnitude trend. 
However, we caution that firm conclusions are difficult 
to draw, since at least part of the scatter is likely the 
result of larger observational errors, given the difficulty 
of measuring accurate photometric parameters for nuclei 
residing in luminous, high surface brightness galaxies. 

Figure [12] also shows the weighted line of best fit for 
nuclei in host galaxies fainter than Bt = 13.5: 

(g - z) DUC = -(0.059 ± 0.034) g DUC + (2.1 ± 0.7). (10) 

Such color-magnitude (or possibly metallictiy-mass) rela- 
tions are generally thought to be a sign o f self-enrichment 
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ailin fc Harris! 120091) . It would not be surprising to 
observe the same self-enrichment in nuclei, given the 
location of the nuclei at the centers of their host galaxies, 
where compressive tidal forces would aid in the retention 
of chemically enriched gas. 

The colors of the nuclei compared to the mean color of 
their host galaxy 's GCs (calculated using the GC sample 
from Paper VIII) are examined in the right-hand panel 
of Figure [13] We find only a very weak trend that 
redder nuclei also have redder GCs, where a weighted 
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Fig. 12. — Color-magnitude diagram for the 31 nuclei identified 
in this study, with colors derived using 4-pixel apertures. Point 
size is scaled with magnitude as indicated. The dashed line is the 
weighted best fit for galaxies fainter than B-r = 13.5. Galaxies 
with B-r < 13.5 or with unusually red or blue nuclei are labeled. 
The mean and standard error of the mean for each luminosity bin 
are indicated by the outlined squares. Right: Histogram of the 
nuclei colors showing a possible bimodal, or skewed, distribution. 



Fig. 14. — Top: Galaxy colors (blue open squares, from lPaper VI) 
and nucleus aperture colors (black closed circles) plotted against 
host galaxy magnitude. Only galaxies that we find to be nucleated 
are shown. The black short-dashed and blue long-dashed lines 
are the best fit to the galaxies and nuclei, respectively. Bottom: 
Difference between galaxy and nucleus color, as a function of host 
galaxy magnitude. The black dotted line marks a difference of zero, 
while the red short-dashed line shows the best fit to all points. On 
average, nuclei are ~ 0.3 mag bluer than their host galaxies. 




fe-z) (gal) (mag) (g-z)(glob) (mag) 

Fig. 13. — Nucleu s aperture colors plotted against host galaxy 
colors fromlPaper VI (left) and mean colors of the GC sample from 
Paper VIlJ (right). The sizes of the circles are proportional to 
the magnitude of the host galaxy. The dotted lines indicate equal 
colors, while the dashed lines show the weighted best fit relation 
for the plotted points. 

least-squares fit gives (with standard errors) 

(g-z)nuc = (0.64±0.84) ( ff -z) g i ob +(0.30±0.83). (11) 

Since mean GC color has been found to correlate 
with that of the host galaxy (e.g., ILarsen et al.l 120011 ; 
IPeng et al.ll2006h . we might expect to find a relation be- 
tween the colors of nuclei and their GCs, given that we 
also find a correlation between nuclei and galaxy colors, 
plotted in the left-hand panel of Figure [13] The weighted 



best-fit line with standard errors is given by 

(g-z) auc = (1.84 ±0.39) (g-z) sal + (1.29 ±0.47), (12) 

which indicates that bluer nuclei tend to lie in bluer host 
galaxies, and vice-versa. The nuclei are also found to 
have a larger range in colors, and are in most cases bluer, 
than their host galaxies. 

In Figure 1 141 we show galaxy and nuclei colors as a 
function of host galaxy luminosity. The colors of both 
the galaxies and the nuclei are found to become redder 
with increasing host luminosity: 

(ff-z)gai = -(0.072 ± 0.010) B T + (2.3 ±0.1), 
(ff-z)„uc = -(0.13 ±0.04)^ + (2.8 ±0.5). 

(13) 

where the errors on the fitted parameters are the stan- 
dard errors. We find the nuclei colors to vary more 
steeply with host luminosity than those of the galaxies, 
although the trend for the nuclei is quite weak for galax- 
ies fainter than Bt ~ 13. Examining the offset between 
galaxy and nucleus colors reveals that those nuclei that 
are redder than their hosts lie predominantly in high- 
luminosity galaxies. The weighted least-squares relation 
and standard errors for the color difference is given by 

A [g _ z) = (0.056 ± 0.033) B T - (0.54 ± 0.48). (14) 

On average, we find the nuclei to be bluer than their 
hosts by (A(g - z)) = 0.28 ± 0.04 mag. If we exclude the 
nuclei in galaxies with Bt < 13 (the regime in which the 
nuclei are found to be redder than their hosts) we obtain 
a mean offset of (A(<? — z)) = 0.32 ± 0.03 mag, where the 
errors are the standard error of the mean. 

4. DISCUSSION 
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4.1. The Role of Environment: Comparison to the 
ACSVCS 

As described in our Fornax survey was preceded by 
a similar study of 100 early-type g alaxies in the Virgo 
cluster (ACSVCS, iCote et al.H2004l) where an investiga- 
tion into the propertie s of t he nuclei in ACSVCS galax- 
ies was carried out by IC06I . Our prime motivation for a 
study of galaxies in the Fornax cluster is to provide a first 
glimpse into the properties of nuclei in two, rather differ- 
ent, clusters, and an assessment of the role played by en- 
vironment in nucleus formation and evol ution. The inter- 
ested reader is referred to § 1 of iPaper 1 . which compares 
some key properties of the two clusters. Briefly, Virgo is 
overall a much larger cluster, with a mass almost 10 times 
that of Fornax (M 2 nn ~ 4.2 x 10 14 vs. (1-7) x 10 13 
(iMcLauehlinl I199& iTonrv et al.l 120001 : iDrinkwater et al.l 
120011) 1. and a velo city dispersion twice as large (a v ~ 760 
vs. 3 74 km s -1 (|Binggeli et al.l 119871 : IDrinkwater et al.l 
I2001D ). Compared to the Virgo Cluster, Fornax is poorer 
(Rich ness Class vs. 1, (jAbell et alj|1989t iGirardi et al.l 
1995)) and more compact (R200 ~ 0.7 vs. 1.55 Mpc). 
Its intracluster medium (ICM) has both lower tempera- 
ture (1 20 vs. 2.58 key) and metallicity (0.23 vs. 0.34 
solar) (jFukazawa et al.|[l998|) . with the Fornax electron 
density at a given radius being about 1/4 that of Virgo 
(jNulsen fc BohringeiifT995HPaolillo et al.ll2002fl . 

In this section, we will direct ly compare the results 
from both surveys. While C06 used King profiles for 
the nuclei in their paper, the ACSVCS results have since 
been updated with Sersic model fits to the nuclei, which 
allows a fair compariso n between the two studies|3 Dis- 
tances from IPaper VI were used to calculate absolute 
magnitudes and physical sizes for both Fornax and Virgo 
galaxies. We note that the two surveys have slightly dif- 
ferent cutoff magnitudes (Bt — 16 for Virgo and ~ 15.5 
for Fornax) and that the distance mo dulus of F ornax is 
^0.5 mag larger than that of Virgo (jPaper VI ). There- 
fore, the Virgo galaxies can reach absolute magnitudes 
roughly 1 mag fainter than those in Fornax. 

4.1.1. Frequency of Nucleation 

In Figure EH we plot the frequency of nucleation of 
the Virgo and Fornax program galaxies as a function of 
their absolute blue magnitude. The Virgo galaxies ap- 
pear in red, and the Fornax galaxies are shown in blue. 
In the upper panel, we overlay histograms for all galax- 
ies (hatched) and nucleated galaxies (solid). This figure 
demonstrates how the Virgo galaxy magnitudes extend 
to ~ 1 mag below those of Fornax, as explained above. 
Our Virgo sample contains 100 galaxies, 67 of which are 
found to be nucleated, so we obtain a total frequency of 
nucleation, /„ = 67 ±8%. This is in excellent agreement 
with the value of /„ = 72 ± 13% found for our full Fornax 
sample. 

The bottom panel shows the frequency of nucleation in 
each luminosity bin. Both clusters exhibit very similar 
distributions with f n = for the bright galaxies, while 
fainter than Mb ~ —19.5, /„ continuously stays above 
~ 70%. Since our Virgo sample has 84 galaxies below 
Mb = —19.5, and our Fornax sample has 35, we find 
the total frequency of nucleation for galaxies fainter than 




M„ (mag) 

Fig. 15. — Same as Figure [4] but using absolute magnitudes 
and showing both ACSVCS and ACSFCS program galaxies (143 
objects in total). Top: Luminosity distribution of the program 
galaxies for Virgo (solid and hatched red histograms) and Fornax 
(solid and hatched blue histograms). The solid histograms show 
the distribution of the 67 Virgo and 31 Fornax galaxies found to be 
nucleated by the ACSVCS and ACSFCS. Bottom: The percentage 
of galaxies found to be nucleated (fn) for Virgo (red squares) and 
Fornax (blue circles). 

M B = -19.5 to be 80 ± 10% for Virgo and 89 ± 16% for 
Fornax. 

Both IC06I and this study have shown that this sharp 
increase in frequency of nucleation compared to previous 
ground based studies (the VCC and FCC) is due m ainly 
to surface brightness selection (see Figure 7 and 8 in lC06l 
and Figure [5] in this work) , which can be attributed to 
the improved resolution and depth offered by the ACS 
imaging. That is, the excellent angular resolution of HST 
has allowed us to uncover previously undetected nuclei 
in both very high surface brightness galaxies, where the 
nuclei are difficult to distinguish from the main body, and 
low luminosity galaxies, in which the nuclei may lie below 
the magnitude limit of the older photographic surveys. 

4.1.2. Nucleus-to- Galaxy Luminosity Ratio 

As in H3.3I and Figure [71 absolute nucleus magnitude 
has been plotted against absolute galaxy magnitude in 
the top panels of Figure UHl Relations of the form 
Equation [7J have been fitted using weighted least-squares 
to the Virgo and Fornax samples, both separately and 
combined, and the parameters and standard errors are 
recorded in Table [SJ the results of which are in agree- 
ment to within the errors for both galaxy samples. 

We also plot nucleus-to-galaxy luminosity ratio r\ as 
a function of absolute galaxy magnitude in the bottom 
panels of Figure [7] The values for the mean and stan- 
dard deviation of r\ are given in Table [5j Taking the mean 
nucleus-to-galaxy luminosity ratio of both data sets com- 
bined, we obtain the following values for each band: 

(rig) = 0.37% ±0.04% , . 

(n z ) = 0.34% ±0.04%, ^1 

which gives a mean value for both bands of 



Sec https : //www. astrosci . ca/users/VCSFCS/Data_Products .html 



(77) = 0.36% ± 0.03%. 



(16) 
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Fig. 16. — Same as Figure[7] but using absolute magnitudes and 
including 67 ACSVCS and 31 ACSFCS nuclei. Top: Nucleus mag- 
nitude plotted against host galaxy magnitude, for the Virgo (red 
squares) and Fornax (blue circles) galaxies found to be nucleated, 
in the g- {left) and z-bands {right). The lines show the weighted 
best fit relations, with the slope held fixed at unity {solid) and 
allowed to vary {dashed). The red and blue lines correspond to fits 
to the Virgo and Fornax samples respectively, while the black lines 
show the fits to the combined sample. Bottom: Nucleus-to-galaxy 
luminosity ratio, n, against host galaxy magnitude, for the g-band 
{left) and z-band {right). The solid and dotted lines show the 
mean and one standard deviation, respectively, while the dashed 
line shows the best fit relation given by the dashed line in the upper 
panel, recast in terms of log(rj) and host magnitude. 
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Virgo and Fornax Nucleus-to-Galaxy Luminosity Ratios 
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TABLE 6 

Virgo and Fornax Nucleus Luminosity Function 



Sample 


Bandpass 


7Tl 

(ma"g) 




(m?g) 


ACSFCS 
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-11.54 ± 0.03 
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58 ± 0.02 


ACSVCS 
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-11.45 ± 0.02 
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87 ± 0.02 


ACSFCS 




-12.67 ± 0.03 


1 


70 ± 0.02 


ACSVCS 




-12.80 ± 0.02 


2 


16 ± 0.02 



The quoted errors refer the standard error on the mean. 

Finally, we note that, due to the definition of 77, the 
best-fit relation from Equation[7]can be recast in terms of 
log(?7) and galaxy magnitude, where a v = —0.4 (ai — 1) 
and ft v = — 0.4/3i. This relation is plotted as the dashed 
line in the bottom panels of Figure [3 and we find that 
we do not see any significant trend between nucleus-to- 
galaxy luminosity ratio and galaxy magnitude. 

4.1.3. Nucleus Luminosities and Sizes 

In Figure [T71 histograms of nuclei luminosities for both 
our Virgo and Fornax sample are compared. The param- 
eters of the weighted maximum-likelihood fit of a normal- 
ized Gaussian to each sample are given in Table [51 where 
the errors on the fitted parameters are the standard er- 
rors. Although we find differences between m° between 
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Fig. 17. — Same as Figure[8] but using absolute magnitudes and 
including 67 ACSVCS and 31 ACSFCS nuclei. The luminosity 
functions for both the Virgo (red squares) and Fornax (blue circles) 
nuclei are shown, in the g-band {top) and z-band {bottom). Both 
data sets have been fitted with a normalized Gaussian. 
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Fig. 18. — Same as Figure ITOl but using parsecs and including 
67 ACSVCS and 31 ACSFCS nuclei. The distribution of half-light 
radii for both Virgo (red squares) and Fornax (blue circles) nuclei 
are shown. The red and blue vertical dotted lines indicate the 
adopted resolution limit of ~ 0'.'025, which corresponds to 2.0 pc 
in the ACSVCS and 2.4 pc in the ACSFCS. 



the two surveys, there amounts are comparable to the 
errors estimated for the nuclei magnitudes. 

We plot a histogram of nucleus sizes in Figure [T5] for 
both our Virgo and Fornax samples. Although there is a 
large range in size (the very large Virgo nucleus belongs 
to VCC 1178), most nuclei appear to have radii < 10 
pc. The typical sizes are in good agreement, with me- 
dian values of 5.7 pc in the g-band for both clusters, and 
7.2 pc and 7.0 pc in the z-band for Virgo and Fornax 
respectively. 

4.1.4. Other Properties 
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In addition to the above properties, we find the Virgo 
and Fornax nuclei to be remarkably similar in a num- 
ber of other ways. First, and most obviously, both the 
ACSFCS and ACSVCS galaxies exhibit a trend along the 
luminosity function in which their central surface bright- 
ness profiles gradually change from having a lumi nosity 
"deficit" to an "excess": see, e.g., Figures 3 and 4 m lCM 
Figure 1 of iPaper U . Figures [2] an d [3] here, a s well as a 
detailed discussion of this trend in IPaper IVl . 

Plotting surface brightness against magnitude, the nu- 
clei are found to hav e diffe rent scaling relations than the 
GCs (see Figure 18 in lC06l and FigurefTTIhere). Although 
IC06I used integrated nucleus colors in their study, our use 
of aperture colors may be a more appropriate compari- 
son to the King profiles used to determine the integrated 
nucleus magnitudes. Indeed, the best-fit line parameters 
outlining the color-magnitude relat ion fo r the nuclei with 
Bt < 13.5, given by Equation 13 in lC06l and Equation ITOl 
in this work, are in good agreement. 

Overall, we find a striking similarity between the nu- 
clei of Virgo and Fornax, despite the clear environmen- 
tal differences between the two clusters. This agreement 
suggests that the physical characteristics of individual 
galaxy clusters (such as ICM density), or the processes 
that depend on them (such as ram pressure stripping ef- 
ficiency), do not play a dominant role in the formation 
and evolution of nuclei in early-type galaxies. Thus, it 
seems we can consider the nuclei examined here as being 
representative of those in early-type galaxies in general. 



4.2. Extension to Low Luminosity: Comparison to the 
Local Group 

Although the ACSVCS and ACSFCS provide a reli- 
able measurement of the nucleation frequency for galax- 
ies brighter than Mb < —15 mag, it is instructive to 
consider the frequency of nucleation in galaxies fainter 
than this magnitude limit. We can do so by examin- 
ing the members of the Local Group, where the small- 
est observed dwarf galaxies reach magnitudes faint as 
My = — 1.5 mag an d can have effective radii on the order 
of ~ 30 pc (see, e.g. JMartin et al.ll2008|) . As sample com- 
pleteness is a concern for such faint, compact systems, we 
focus on the subset of early-type galaxies brighter than 
Mb ~ -8. 

At present, there are 25 known early-type galaxies in 
the Local Gro up brighter than thi s limit (compiled from 
Mated 119981: IMcConnachie et al.l 120051: iBrasseur et""aH 
20111) . Of these, only two (NGC 205 and M32) 



are brighter than t he ACSVCS limiting magnitude of 
M B = -15 mag ()Mated Il9£)8l). both o f which are 
known to be nucleated ( e.g. iKentl 119871 : lLauer et al.1 



19981: iMated 119981: iButler fc Martmez-Delgadd I200E 



De Riicke et al. 2006). Moving down the luminos- 
ity function, at most six other galaxies may con- 
tain either nuclei or kinematically /structurally dis- 
tinct features near their core, listed in order of de - 
creasing lum inosity: NGC 147 (iDe Riicke et al.l f2006). 
Sagittarius (IMated 119981: ILavden fc Saraiedinil 1200(1 



iMonaco et all 



dColeman et al.l 12004. 120051: iColeman fc de Jong] l2008h. 
Sextans dKlevna et al.l 12004 iWalker et all 12009). An- 
dromeda II ()McConnachie &: Irwin! 120061 ). and Ursa Mi- 



200."). jBellazzini et al.l 120081), Forna x 



nor (|Klevna et al.|[200l iPalina et al.|[200i r ro l Consider- 
able caution is advisable here since, in some cases (e.g., 
in Ursa Minor and, especially, in Sextans), the "nuclei" 
are rather subtle substructures (sometimes only appar- 
ent with the addition of kinematic data) that bear little 
resemblance to the prominent, compact nuclei seen in the 
faintest ACSVCS and ACSFCS galaxies. Yet, even with 
this liberal definition of a "nucleus", only 8 out of the 
27 Local Group early-type galaxies (/„ = 30%) can be 
classified as nucleated. If we exclude Fornax, Ursa Mi- 
nor and Sextans from the list of nucleated galaxies, then 
/„ falls to 19%. While it is possible that some nuclei 
have yet to be discovered, it seems certain that many of 
the f aint Local Group galax i es do not contain a nucleus; 
e.g., IMcConnachie &: Irwin! (|2006l ) studied of structural 
properties of six Andromeda satellites using deep, homo- 
geneous imaging, and found a nucleus in only a single 
object (And II). 

We conclude that the frequency of nucleation along 
the Local Group sample is clearly far lower than in ei- 
ther our ACS surveys of the Fornax or Virgo clusters. 
Why is there such a large disparity in /„? We specu- 
late that the lack of nuclei in very faint galaxies could 
be related to the general absence of GCs in galaxies be- 
low M B 12 (see, e.g.. iPeng et~aTll2008D . If nuclei in 

low-mass galaxies are assembled through GC infall and 
mergers (see ^4.3.11) . then the presence of GCs would 
obviously be a prerequisite for nucleus formation. The 
faintest galaxies in the Local Group known to contain 
GCs are Sagittariu s and Fornax , with Mb = —12.8 and 
— 12.6 respectively (Mateo 1998). The former is unques- 
tionably nucleated, while Fornax may meet the defini- 
tion of a nucleated galaxy (see above). Because no Local 
Group dwarfs below this magnitude are known to con- 
tain GCs, such galaxies might have been unable to form 
a nucleus if star cluster infall is the dominant mode of 
nucleus formation in low-mass systems. 

It is also interesting to note that, assuming a constant 
nucleus-to-galaxy luminosity ratio of 0.4%, then the ex- 
pected nucleus magnitude of a Mb = —12.6 mag host 
would be Mb = —6.6 mag. This corresponds closely to 
the mean turnover magnitude of the globu lar cluster lu- 
mino s ity function My sa —7.5 mag (e.g., Jacobv et aLI 
1992; Harris 200fl IBrodie fc Straderl 12006). suggesting 
that galaxies may be unable to form nuclei at the point 
where the expected nucleus luminosity would fall below 
the typical GC luminosity. 

However, as caveats we firstly n ote that the nucleus 
of Sagittarius (M onaco et al.l 120091) as well as the very 
central region of the Andromeda satellite NGC 205 
(|Siegel et al.l 120071) have been observed to have under- 
gone multiple star formation episodes, which indicates 
that other processes in addition to GC accretion must 
have shaped their formation history. In addition, the nu- 
clei late-type dwarfs have been shown to not fo rm form 
exclusively fr om GC infall (iWalcher et al.ll2006l e.g.,) or 
gas accretion (|Hartmann et al.ll201lD . even though it has 
been observed that GC specific frequency is independent 

19 Although it is traditionally classified as non-nucleated, we 
include the Fornax dSph in this list since GC #4 is located ~ 
half a core radius from t he galaxy photocenter (see Figure 1 of 
Coleman & de Jong 2008) and might thus be classified as a dwarf 
with an offset nucleus if moved to the distance of the Virgo or 
Fornax clusters. 
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of morphology (|Georgiev et a l. 2010) and thus should be 
the same for both early- and late-type dwarfs. 

4.3. Formation and Evolution Models 

The origin of nuclei remains an open theoretical 
problem, with two main avenues of nucleus formation 
presently considered most viable. The first proposes that 
a galaxy's star clusters will experience orbital decay due 
to dynamical friction and spiral inwards, eventually coa- 
lescing at the center of the galaxy. The second formation 
mode focuses on gas accretion at the center of the galaxy, 
followed by star formation. Some similarities in the scal- 
ing relations of nuclei and black holes (see §2} have also 
given rise to models that consider the formation of both 
types objects in a shared context. In this section, we 
shall examine theoretical studies of nuclei formation in 
light of our new results, as well as models that explore 
the relationship between nuclei and black holes. 

4.3.1. Disstpationless Infall of Star Clusters 

iTremaine et al.l (|1975f) first suggested that the nucleus 
of M31 was formed from GCs that spiraled inward to the 
galaxy center due to dynamical friction, and this mecha- 
nism continues to offer an attractive explanation for the 
assembly of nuclei in at least some galaxies. Of course, 
not all clusters that come close to the center of a galaxy 
will necessarily contri bute to the format i on, or growth, 
of a stellar nucleus; as ICapuzzo-Dolcettal (|1993f ) showed, 
dynamical friction and tidal stripping are competitive 
processes, where GCs are more readily destroyed by large 
nuclei, limiting nucleus growth. 

Nevertheless, some fraction of GCs are expected to 
avoid tidal disruption and could contribute to cither 
nucleus formation, or the growth of pre-existing nu- 
cl ei. Evidence in favor of th i s proc ess was described 
in iCapuzzo-Dolcetta fc Tesserl (|1999ft , who pointed out 
that the radial distribution of GCs in galaxies is 
less centrally concentr a ted t h an the halo sta r s (see 
also IMcLaughlinl [19951 Il999l : ICote et al.l l200ll [20031 : 
IPeng et al.l 120081) . Such "missing" clusters could have 
contributed to nucleus formatio n. Monte Carlo simu- 
lations based on this premise by lLotz et all (|2001l) pre- 
dicted nuclei luminosities for dEs with — 17 < Mb < — 12 
that were consistent with observations for the brighter 
galaxies within this range, although they were overes- 
timated for less luminous ones. The over-prediction of 
nuclear luminosities in their low-mass systems resulted 
from their short dynamical times — meaning that nuclei 
are able to grow very efficiently — in spite of the fact 
that these galaxies ha ve relatively few star clusters (see 
e.g., IPeng et alJl2008l ). 

Numerical simulations by lOh fe Linl (|2000f l and 
similar, higher resolution N-bod y simulations by 
ICapuzzo-Dolcetta fe Miocchil (|2008allbT ) were able to suc- 
cessfully reproduce the observed surface brightness pro- 
files of known nucleated galaxi es. A dependenc e on lo- 
cal tidal field was found in the IQh fc Linl (pOOOh model, 
where disruptive tidal forces on the outskirts of galaxy 
clusters would alter GC orbits, increasing dynamical fric- 
tion timescales and decreasing nucleation frequency. The 
ICapuzzo-Dolcetta fc Miocchil (|2008al fbh models suggest 
that, if linear scaling is assumed, then the observed nuclei 
could have formed from the infall of tens, to hundreds, 
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Fig. 19. — Top: Histogram of masses for the 143 galaxies from 
the ACSVCS and ACSFCS surveys. Bottom: Dynamical friction 
timescales, Tqf, plotted as a function of galaxy mass. Two sets 
of curves are shown. The dashed blue curves show calculations 
for initial GC radii, Rj , eq ual to the galaxy effective radii (see the 
lower panel of Figure 1116 , while the dotted red curve shows Ri 
fixed to 1.3 kpc, the median effective radius for ACSFCS galaxies. 
In both cases, Tdf is plotted for five GC masses: 0.1,0.25, 0.5, 1.0 
and 2.0 million solar masses. Note the sharp decline in Tdf for 
low-mass galaxies. 

of GCs (see also §4.9 and §5.2.4 of [C06h. Both simu- 
lations found that nuclei may begin to coalesce away 
from the galaxy photoccntcr, although to quite differ- 
ent extents: i.e., up to ~ 0.3 kpc and settling within 
~ 1 Gyr in IQh fc Linl j OOO) , and ~ 4 pc away in 
ICapuzzo-Dolcetta fc Miocchil (l2008a|lbT> . 

Other simulations bv lBekki et al.l ()2004f ) observed that 
the scaling relations of nuclei formed through mergers of 
GCs would be notably different than those of the GCs. 
In $3.51 we discussed that the predicted scaling relation 
for nuclei in these simulations, R e oc was gener- 

ally in good agreement with observa tions (see Figure ITT|). 
More recent work by [Bckki (20l3) focused on simula- 
tions of star cluster infall due to dynamical friction in 
disk galaxies. He found that the effectiveness of dynam- 
ical friction did not depend strongly on bulge mass, but 
increased with smaller disk mass, and with larger disk 
mass fraction, galaxy surface brightness, and star cluster 
mass. The ratio of nucleus mass to disk mass was found 
to decrease as a function of increasing disk mass, with 
a mass ratio of > 0.4% for smaller disks, and < 0.1% 
for disks with masses M > 10 9 Af Q . However, star clus- 
ter mergers on to a disk may not be sufficient to explain 
nuclei fo rmed in My ~ —19.5 spirals. N-body simula- 
tions bv lHartmann et all (|2011h . which aimed to repro- 
duce the observed kinematics of the nuclei in M33 and 
NGC 4244, found that star cluster accretion on to a disk 
did not produce the necessary line-of-sight velocity rise, 
and at least half of the nucleus mass had to come from 
gas dissipation. 

Some provisional evidence for dissipationl ess formation 
in at least some galaxies was presented in iPaudel et al.l 
(|2011[ ). who used optical spectroscopy for Virgo cluster 
dwarfs to study both their stellar populations and those 
of their nuclei. Despite the small sample and the different 
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environment (Virgo vs. Fornax), their data present an 
interesting opportunity to speculate on p ossible forma- 
tion m echanisms for the ACSFCS nuclei. iPaudel et al.l 
(j 2 1 lh found that nuclei in a handful (5) of the faint 
(—16 < Mb < —14) galaxies in their sample were older 
and more metal poor than their hosts, which is certainly 
suggestive of a connection to GCs. At higher luminosi- 
ties, most of their nuclei were found to be younger than 
their hosts. While inconsistent with nucleus formation 
from old globular clusters, this observation may still be 
compatible with cluster infall, as our observations and 
many others have shown that ongoing star cluster for- 
mation can be present throughout so me galaxies (e.g., 
I Anders et al1l200i iKveong et al.ll2010h . In the ACSFCS 
sample, FCC 119, FCC 90 and FCC 26 are possible ex- 
amples of Mb > —19.5 galaxies with young cluster sys- 
tems. 

Additional support for such a scenario may come from 
the GC luminosity functions in Virgo and Fornax galax- 
ies. The widths of GC luminosity functions are known 
to decrease significantly with galaxy luminosity, a trend 
that is accompanied by a sl i ght decrease in turnover mass 
(| Jordan et al.ll2006l [20075 (paper Villi ). This truncation 
of the GC population on the bright end of the luminosity 
function may be caused, at least in part, by the shorter 
dynamical friction times as galaxies become less massive, 
although ot her (external) proces ses could also play a role 
(see §7.2 of I Jordan et"aLll2007bt ). 

We revisit the question of star cluster infall efficiency 
by calculating the dynamical friction timescale, Tdf, for 
all galaxies in our ACS surveys of Fornax and Virgo. 
The upper panel of Figure [TO] shows the distribution 
of galaxy masses from the combined surveys (filled his- 
togram), while the lower panel show the dependence of 
Tdf on galaxy mass, which is given by 



Tdf = 



2.64 x 10 2 



In A 



2 kpc 



250 km s~ 



10 6 M ( 



GC ' 



Gyr 



(16) 

Here Ri is the initial galactocentric radius of the star 
cluster, v c is the circular velocity of the (assumed isother- 
mal) galaxy, and Mac is the mass of the star cluster 
(|Binnev fc Tremaind 120081 ). In this equation, In A is the 
coulomb logarithm, which is defined as 
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where & max is the maximum impact parameter between 
the cluster and the interacting particle (a star of mass 
m). Following lLotz et al.l (j200l|), we assume v c ~ \f2a 
where a is the integrated-light velocity dispersion mea- 
sured within R e / A from McLaughlin et al. (2012, in 
prep.). We also take 6 max = R e for all galaxies, with 
R e measured directly from the ACS imaging (see WS.5\ 
and Figure 1TT1) . 

Calculations have been carried out for five different star 
cluster masses (i.e., 0.1, 0.25, 0.5, 1 and 2 million solar 
masses and for two assumpti ons for Rj. In t he first 
case, we take Ri — R e (see also lLotz et al.|[200ll ) which 

20 Recall that in the Milky Way, the GC mass corr esponding to 
the p eak of the luminosity function is 2.4 X 10 5 Mq (McLaughlin 
IT999D . 



is shown as the dashed blue curves in Figure [H)l In the 
second case, we simply fix Ri at the median effective ra- 
dius, 1.3 kpc, for all galaxies in the ACSFCS sample. 
The results in this case are indicated by the dotted red 
curves in Figure [19] Although Tdf clearly varies with 
the assumed cluster mass and the precise choice of Ri, 
the strong mass dependence noted by previous investiga- 
tors is clearly apparent in this figure. In particular, the 
dynamical friction timescales are dramatically shorter in 
galaxies with Ai* < 10 w M.q compared to higher-mass 
galaxies. We conclude that star cluster infall seems like a 
viable, indeed a likely, candidate for the growth of nuclei 
in low- and intermediate-mass galaxies in our sample. 
For the highest-mass galaxies, the mechanism appears 
much less viable given the fact that, in these systems, 
Tdf greatly exceeds the Hubble Time for all but the 
most massive and centrally concentrated star clusters. 

Finally, we conclude this section with some final re- 
marks on Figure 1111 which compared the structural pa- 
rameters of nuclei to those of GCs and their host galaxies. 
While there is, as noted in £ 13.51 good agr eement with 
the nu clei size-mass relationship found by iBekki et al.l 
(2004) from simulations of GC mergers, there are rea- 
sons to believe that a single relation cannot be appropri- 
ate for all nuclei which, in our sample, span more than 
four de cades in mas s . For comparison, the simulated nu- 
clei of IBekki et all (l200l span a factor of just ten in 
mass. It is to be expected that the precise form of the 
size-mass relation in the context of the GC merger model 
will be different in different mass regimes. For instance, 
when only a small number of mergers contribute to the 
nucleus, we expect from the virial theorem and conser- 
vation of energy that R e cx M®' 5 . At later times, when 
the mass of the nucleus greatly exceeds the mass of an 
accreted GC, the relation should steepen to R e oc Al*. 
These scaling relations, shown in the lower panel of Fig- 
ure 1111 are in good agreement with the observed sizes 
and masses. 

All in all, based on the existing data, we believe that 
cluster infall must have played an important role in the 
formation of the nuclei the low- and intermediate-mass 
hosts within our sample. At the same time, the red col- 
ors of some of the largest and most massive nuclei ( £13.61) 
present a strong challenge to this model, suggesting that 
an additional process — most likely the dissipational in- 
fall of metal-rich gas — likely begins to dominate the for- 
mation of nuclei in galaxies of progressively larger masses 
dMihqs fc Hernquistl[l99l ICOa iPaper ID : iHopkins et al.l 
2i HIS. 120091 ) ~ 

4.3.2. Dissipational Infall of Gas 

It has long been suspected that nuclei could form 
through star fo rmation following th e accretion of gas in 
galaxy centers (jvan den Bcrgh 1986), although the exact 
origin of the gas, and the mechanism that triggers the 
inflow, are matters of debate. 

In some models, the gas is assumed to orig inate from 
outside the galaxy. iDavies &: Phillippsl (|1988| ) proposed 
that dEs may be formed from fading stellar populations 
in dwarf irregulars, where the accretion of Hi gas induced 
starbursts, the fin al one occurring in the center and form- 
ing the nucleus. I Silk et all (jl987) predicted that the in- 
tergalactic medium (IGM) could fall into dwarf galaxies 
when it is cooled and compressed during group forma- 
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tion. This model noted that dwarfs closer to large galax- 
ies may not be able to form nuclei as efficiently, since the 
large galaxy's ti dal field makes it diffic ult for the dwarf to 
capture the gas. iBabul fc Reesl (|1992f l found an opposite 
trend with environment: they observe that nucleus evo- 
lution may depend on local IGM density, because this 
determines whether supernova-driven gas outflows are 
able to escape. Dwarfs in low-pressure regions would 
have their gas ejected and then fade away, while winds in 
dense environments would be restricted to the starburst 
region by the IGM. This confinement could cause gas to 
cool and recollapse, creating two short or one prolonged 
starburst. 

Gas might also be funneled to the centers of 
galaxies which have disks and axisymmetric features. 
IMilosavlievicI (|2004D suggested that in spiral galaxies, 
magneto-rotat ional instabi l ity in the disk transport s gas 
to the center. IBekki et all J2006) and IBekkil (|2Q07f ) per- 
formed chemodynamic simulations of the inner 1 kpc of 
dwarf galaxies with stellar masses of 2.5 7 < M sp h < 
1.0 x 10 9 , to explore the remnant created through dis- 
sipative merging of stellar and gaseous clumps formed 
from nuclear gaseous spiral arms in a gas disk. The sim- 
ulations produced nuclei which that rotating and flat- 
tened, consisting of stars with varying ages and metal- 
licities. Although the initial clump was found to form 
off-center (about 200 pc by visual inspection of the simu- 
lation data), it would fall into the center within 100 Myr. 
They found that overall, the nuclei were characteristi- 
cally younger and more metal rich than the host, with 
more massive hosts creating more metal-rich nuclei. Gas 
settling timescales increased with decreasing dwarf mass 
(due to feedback being more effective in smaller galax- 
ies), so low mass dwarfs were found to have younger and 
bluer nuclei. More massive and dense nuclei were formed 
in more massive dwarfs with deeper central potentials, 
and both the mass and mass fraction of the nucleus were 
found to increase with spheroid mass. Nuclei in high sur- 
face brightness galaxies should also have higher surface 
brightness, owing to the increased dynamical friction due 
to higher stellar densities. The nucleus surface brightness 
was strongly dependent on the gas fraction of the host, 
and thus may be more likely to form in this manner in 
late- type galaxies with relatively large amounts of gas. 
Finally, the addition of a central black hole to the simu- 
lation had little effect on the properties of the remnant 
nucleus. 

Anoth er source of nu clear material, which was first pro- 
posed bv lBailevi (|1980D . could arise from stellar winds. It 
was found that only a small (~ 1O 6 M ) amount of gas 
was needed to cause an inflow fo r an elliptical galaxy 
with M ga i - 10 n M Q . iSethl (|2010[ ) observed that such a 
mechanism could produce the age, abundance gradient, 
and rotation curve seen in the nucleus of M32. 

The dissipative infall of gas to th e galaxy center can 
also b e induced by galaxy mergers. Mihos fc Hernquistl 
(1994) performed N-body simulations of disk galaxy 
mergers, where they found that gas dissipation and 
the star formation that followed created dense stellar 
cores in the remnant. S imilar highe r resolution simula- 
tions bv Hopkins et al.l (|2008l I2009D . showed that grav- 
itational torques during gas-rich mergers removed the 
angular momentum of the gas, which would then un- 



dergo gravitational collapse. The amount of gas in- 
fall was found to largely depend on the progenitor 
galaxy gas fraction, while the addition of a central 
black hole was not found to have a significant effect 
on the the properties of the final remnant. Unfortu- 
nately, these models lacked the resolution to study typ- 
ical nuclei, particularly those in the low-mass galax- 
ies: i.e., apart from a small number of cE galaxies in 
the ACSVCS sampl e, which have likely been heavily 
tidally stripped ( e .g.. lFaberll!973l: IFerrarese et alJfe OOBb: 
| C6te et all 120081 : iChilingarian et al l 120091: iHuxor et al.l 
120111 : McL aughlin et al. 2012, in prep.), the simulated 
galaxies of (|Hopkins et a l. 2009) have masses > 10 10 M & , 
more than ten times larger than the masses of the faintest 
galaxies in the ACS surveys. However, in this restricted 
mass regime, the properties of these simulated galaxies 
are in good agreement with our ACSFCS (and ACSVCS) 
observations. 

Likewise, t he simulations of IBekki et"aLl ()2006D and 
IBekkil (|2007f ). which instead focused on the low-mass 
galaxies, also appear to be consistent with observations, 
including those from our HST/ACS imaging and re- 
sults from ground-based spectroscopy. First, the nu- 
clei in these simulations were found to be younger and 
more metal rich than their hosts, with nucleus metal- 
lic ity increasing with host mass, a trend that was seen 
in IPaudel et al.l f|201 If ) . Second, their finding that low- 
mass dwarfs have younger and bluer nuclei is consis tent 
with some of the nuclei from IPaudel et al.l (|2011l ). as 
well as with the nucleus colors observed in our study. 
Finally, they also found that the mass fraction of the 
nucleus increased with host spheroid mass, and that 
their simulated surface brightness profiles showed nu- 
clei which become more prominent with increasing dwarf 
mass, whereas in low-mass dwarfs the nuclei were barely 
distinguishable. It is therefore possible, as discussed in 
^4.3. 11 that nucleus formation through gas infall may be 
most significant for intermediate- and high-m ass galax- 
ies. In their analysis of the ACSVCS, IC06I noted that 
some of the reddest and brightest nuclei "may be candi- 
dates for th e dense stellar cores that form in numerical 
simulations (|Mihos fc Hernq uist 1994]) when (chemically 
enriched) gas is driven inward, perhaps as a result of 
mergers." Such a result can be reconciled with our nearly 
constant nucleus-to-galaxy luminosity ratio if star clus- 
ter infall accounts mainly for nucleus build up in lower- 
mass galaxies. At intermediate masses, both processes 
could contribute significantly to the growth of nuclei; 
candidates for such hybrid nuclei in the ACSFCS include 
FCC 43, FCC 249, FCC 310, FCC 148 and FCC 301, 
which may consist of both compact and extended com- 
ponents. 

This basic scenario is also consistent with the gen- 
eral view that mergers (which can drive gas to the 
central regions of a galaxy) become increasingly im- 
portant as galaxy luminosity increases, a consequence 
of the hierarchical merging paradigm. The observa- 
tion that galaxy concentration — parameterized by 
Sersic inde x n — varies smoothly with galaxy lumi- 
nosity (e.g-Uerien & Binggelilll997tlGraham et"aLll2003l; 
IFerrarese et al.ll2006bl : McLaughlin et al. 2012, in prep.; 
see also P2. 1() provides strong supporting evidence for 
this picture, as violent relaxation of merger remnants 
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Fig. 20. — Top: Nuclei luminosity fraction plotted against Sersic 
index of the host galaxy, n ga i. Bottom: Nuclei effective radius as a 
function of n ga j . The dashed line in each panel shows the weighted 
best-fit linear relation; unweighted fits are shown by the dotted 
lines. The nuclei in both clusters show weak trends with Sersic 
index (or, equivalently, galaxy mass) in the sense that the central 
"excess" above the fitted Sersic model seem to be brightest and 
largest in galaxies with the largest n ga [. These galaxies have likely 
undergone fewer mergers and accretions than those with n ga j ~ 1. 



is thought to be respo nsible for the creat ion of de Vau- 
couleurs profiles (e.g. iBarnesI fl988l I1992D . while Sersic 
index of both bulge and disks of spirals has been show n 
to increase after satellite infall (fElichc-Mo ral et aLll2005l ). 
Figure [20] shows the dependence of two fundamental pa- 
rameters for nuclei — luminosity fraction and effective 
radius — against host galaxy Sersic index (Paper III). 
Those galaxies whose internal structure has likely been 
transformed most extensively through mergers, accre- 
tions and harassment (i.e., those galaxies with high Sersic 
indices) tend to have the most luminous and spatially ex- 
tended nuclei (although the trend between n and rj is sta- 
tistically significant only when an unweighted fit is used). 
These trends are generally consistent with an increasing 
importance for gas dissipation as ones moves to higher 
and higher mass galaxies. Nuclei formed through merger- 
driven gas inflow could also be expected to follow a mass- 
radius scaling relation, as iHopkins fc Quataertl (2010) 
found that stellar systems may have a maximum stellar 
surface density, due to feedback from massive stars. 

One complication with the gas inflow model is that 
it obviously requires the presence of gas, which is 
not consistent with the "classical" picture of early- 
type galaxies. However, both low-mass Es and high- 
mass "dEs" are now recognized to be quite com- 
plex, having been found to co ntain dust, spiral 
arms, embedded disks, and bars dJerien et al.l |2000|; 
Barazza et al.ll2002t IDe Riicke et al.ll2003at iLisker et al.l 



2006bt iFerrarese et al.ll2006bf l. as well as counter rotat- 



ing and kinematically decoupled cores (IDe Riicke et al.l 
l2004HThomas et al.ll2061 iChilingarian et al.ll2008fl. and 
ongoing star form a tion (e.g., IDe Riicke et al.l 12 003b; 
Lisker et all l2006al iCote et al.l 120061: iMichielsen et al l 



20071) . These features suggest that a non-negligible frac- 



tion of intermediate-mass galaxies classified as "early" 



types have experienced some level of morphological 
transformation, likely through mergers, a ccretions, or 
interactions with the cluster environment ([Moore et al.l 
ll996HKazantzidis et aLll2011h . 

It is, in fact, possible that the nuclei in some of our 
early-type galaxi es formed in late-type progenit ors. A 
recent finding by Ems ellem fe van de Venl (|2008l ) noted 
that galaxies with Sersic indices of n < 3.5 have com- 
pressive tidal forces in their central regions, with the 
size of the compressive region increasing with decreasing 
Sersic index. Assuming a constant Sersic index of n = 1 , 
the amplitude of the tidal forces was found to scale lin- 
early with galaxy mass, and form a central massive object 
(CMO) with a constant host mass fraction of ~ 0.5%. A 
CMO growing through gas accretion in this way would 
eventually reach a critical density and luminosity, alter- 
ing the galaxy profile such that it no longer has cen- 
tral compressive forces. Compar ison of this theoretical 
threshold nucleus luminosity with lC06l reveals that many 
observed nuclei are much more luminous than would be 
predicted by this model, which suggests that the nuclei in 
early-type galaxies may have formed in some low Sersic 
index, gas-rich progenitors that have since evolved mor- 
phologically. 

High-resolution observations of molecular and neutral 
hydrogen in these galaxies may be able to constrain the 
role of gas inflow and enrichment in nucleus formation, 
since H2 will highlight regions of star formation, while 
HI is a tracer of processes affected by the ICM and grav- 
itational interactions. Subarcsecond-resolution maps of 
molecular starburst gas — by using ALMA to observe the 
CO transitions and EVLA to detect HI through 1.4 GHz 
emission — would allow the relationship between galaxy 
nuclei and molecular gas to be examined in much greater 
detail than is currently possible. 

4.3.3. Possible Connections to Black Holes 

As discussed in § [H recent observations have uncovered 
the coexistence of nuclei and black holes in intermediate- 
mass galaxies, which may have implications for the evo- 
lution of the central regio ns of galaxies. For instance, 
lHopkins~fc Quataert] (j2010f ) performed simulations of gas 
accretion on to a black hole, which they find can form 
a lopsided, eccentric nuclear disk that exerts a strong 
torque on and drives in the remainder of the gas, pro- 
ducing a system much li ke that found in M31. 

Another simulation bv lBekki fc Graham! ([20101 ) exam- 
ined the merging of two nuclei containing black holes, 
and found the dynamical heating of the cluster from the 
black hole binary expelled stars from the center, with 
the final stellar density of the remnant decreased with 
increasing black hole mass fraction. This type of merger 
could produce observed "core" galaxies with larger b lack 
hole s (as originally noted by lEbisuzaki et aTl 119911 : see 
also iMilosavlievic fc Merrittll200ll ). and shape the inner 
regions of intermediate-luminosity galaxies in which a 
nucleus is difficult to distinguish observationally f£ )3.ip . 
Their simulations further showed that if only one nucleus 
had a black hole, the decrease in stellar density of the nu- 
cleus was less pronounced, as most of the heating comes 
from the black hole binary. In mergers where neither nu- 
cleus had a black hole, the stellar density of the nucleus 
increased. 

If black holes do become an increasingly dominant 
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component of the CMO mass budget in high- and 
intermediate-luminosity galaxies, then they could either 
hinder nucleus growth, or lower the density of the nu- 
cleus through mergers until it is destroyed by black hole 
binary feedback. These effects could create the trends in 
intermediate-mass galaxy surface brightness profiles ob- 
served in this study, where the galaxies undergo a tran- 
sition from central light "excess es" to "de ficits" as they 
become more luminous (see also lPaper IVT ). 

5. SUMMARY 

This HST study examined 43 early-type galaxies in the 
Fornax cluster, imaged in the ACS F475W and F850LP 
bands. Our analysis — performed in both one- and two- 
dimensions — extracted photometric and structural pa- 
rameters for 31 compact stellar nuclei in these early- type 
galaxies. The main results are summarized as follows: 

1. We have compared our ID results to those obtained 
by using 2D image modeling techniques, and found 
the extracted nucleus structural parameters to be 
in agreement for both methods. Although 2D fit- 
ting potentially allows for full structural decompo- 
sition of a galaxy, ID methods enable characteriza- 
tion the outer regions with a single surface bright- 
ness profile. We conclude that ID fits are more 
appropriate for our study, since they allow us to 
easily compare nucleus and galaxy parameters in 
an objective and homogeneous way. 

2. We find that 72 ± 13% of the 43 galaxies in our 
sample are nucleated, which is a significant increase 
from ground-based studies. The nuclei — defined 
as a central excess relativ e to t he inward extrap- 
olation of a Scrsic model (C06) — are found ex- 
clusively in galaxies with Mb > —19.5 < 
1O 1O ' 6 A^0), and the frequency of nucleation for 
galaxies fainter than this magnitude is 89 ± 16% 
(31/35). As was found previously in the Virgo 
cluster, nuclei are exceedingly common in low- 
mass, early- type galaxies in the Fornax cluster (i.e., 
M, > 10 9 M Q ). 

3. Most nuclei are not significantly offset from their 
host photocenter — only three are offset by more 
than 0'.'5. We do not find any trend between the 
magnitude of the offset and host galaxy luminosity. 

4. We find a nearly constant nucleus-to-galaxy lumi- 
nosity ratio of ~0.4%. The observed nucleus lu- 
minosity function can be understood therefore in 
terms of the galaxy selection function (and the 
fact that galaxies brighter than Mg < —19.5 do 
not contain nuclei). If we parameterize the nu- 
cleus luminosity function as a normalized Gaus- 
sian, we find peaks at (M g ) = —11.5 and (M z ) = 
— 12.7 mag, which is ~ 40 times more luminous 
than the peak of the GC luminosity function. The 
nuclei are also found to have larger sizes and dif- 
ferent effective surface brightness scaling relations 
than the GCs. 

5. The colors of the nuclei in hosts with Bt < 13.5 
are found to correlate with galaxy colors, as well as 



with galaxy and nucleus luminosities. In particu- 
lar, both the galaxies and the nuclei were observed 
to become increasingly red with increasing galaxy 
luminosity, with the trend being steeper for the nu- 
clei. This leads to a relation between nucleus-and- 
host color difference and host magnitude, where 
nuclei that are more red than their hosts are found 
predominantly in brighter galaxies, and vice versa. 
However, on average most of the nuclei are signif- 
icantly bluer in (g — z) color than their hosts by 
0.28 ± 0.04 mag. 

6. A comparison to lC06l which examined the nuclei of 
early-type galaxies in Virgo, reveals many similar- 
ities between the nuclei in the two environments. 
Both studies find similar frequencies of nucleation 
(increasing sharply from to > 70% for galax- 
ies with Mb > —19.5 mag), surface brightness 
selection effects, nucleus-to-galaxy luminosity ra- 
tios, nucleus luminosity functions, sizes, and color- 
magnitude relations. The trend along the lumi- 
nosity function where the galaxy central surface 
brightness profiles gradually change from having 
a luminosity "deficit" to an "excess" is sh ared by 
both samples (see also lPaper II] ; lPaper IVD . which 
suggests that generic formation and evolution pro- 
cesses largely independent of the galaxy environ- 
ment are involved in shaping the central regions 
of galaxies. Rather, nucleus creation may be more 
contingent on local factors, especially host galaxy 
mass. 

Our conclusion is that, in low-mass galaxies, the dom- 
inant mechanism for nucleus growth is probably infall of 
star clusters through dynamical friction, while at higher 
masses, gas accretion resulting from mergers and torques 
becomes dominant. There is no reason to expect either 
of these processes to be discontinuous, and we argue that 
the relative importance of these processes vary smoothly 
as a function of galaxy mass. We examine the efficiency 
of dynamical friction in our sample galaxies and confirm 
the finding of many previous studies that star cluster in- 
fall is most effective in low-mass galaxies. Based on sim- 
ulations carried out by other researchers, we argue that 
gas infall, followed by central star formation, becomes in- 
creasingly important in high-mass galaxies having Sersic 
indices that may have been inflated by successive merg- 
ers and accretions. There is also some evidence for "hy- 
brid nuclei" in some of the intermediate-mass galaxies in 
our sample: i.e., nuclear components with complex inner 
structures. Simulations that take into account multiple 
formation mechanisms — star cluster infall, gas accretion 
driven by tidal torques and/or accretions and mergers, 
the influence of central black holes, etc — are urgently 
needed to elucidate the processes that drive nucleus for- 
mation in different mass regimes. 

Both dissipationless cluster infall and gas accretion 
models make predictions that nucleus formation would 
depen d on local density (|Oh &: Linl 120001 : iBabul fc Reesl 
Il992| ). Although the fact that we do not find any major 
differences between the nuclei of Virgo and Fornax sug- 
gests that local density may not be a dominant factor 
in their formation, observations that examine the entire 
volume of a galaxy cluster (and that have the sensitivity 
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necessary to detect the nuclei) may help determine the 
role environment plays in shaping the nuclei and their 
hosts. In this context, the forthcoming Next Genera- 
tion Virgo Cluster Survey (Ferrarese et al. 2012), which 
is imaging the entire Virgo cluster to a (10cr) depth of 
g 25.7, should provide important new constraints on 
formation models. 
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APPENDIX 

A. COMPARISON WITH 2-DIMENSIONAL 
SURFACE BRIGHTNESS PROFILE 
FITTING 

<J2] describes the determination of nuclei parameters 
through ID fitting of surface brightness profiles from 
ellipse. This appendix examines the pros and cons 
of ID and 2D methods when measuring parameters for 
the photometric and structural parameters of nuclei in 
early-type galaxies. 

In general, the decision to use a ID or 2D approach 
depends on the specific scientific goals. If a galaxy has 
multiple components (which even for early-type galax- 
ies can include, e.g., bulges, large-scale disks, embedded 
disks, outer/inner rings, bars, shells, dust filaments, dusk 
disks, and even faint spiral armsj3, then, by using 2D 
decomposition, individual structure can, in principle, be 
fitted with separate profiles and the galaxy's composi- 
tion examined in deta i l. Th e 2D fitting program GAL- 
FIT (jPeng et al J 120021 l2010h allows the implementation 
of many surface brightness profile modifications, such as 
variability of their diskiness/boxiness, or the addition of 

21 Although early- type galaxies are often considered structurally 
simple systems, all of these features are fou nd in the sample of 143 
early-type galaxies studied in the ACSVCS (Ferrarese et al. 2006b) 
and ACSFCS (Paper III). 




Fig. 21. — Upper row: F475W image for FCC 190 displayed 
at three different intensity stretches (A, B and C) and two dif- 
ferent magnifications (A/B vs. C). Note the prominent nucleus 
visible in panel (A), and the dramatic changes in ellipticity and 
position angle with radius. Panels (D) and (E): Galaxy model 
constructed using ellipse, with contours overlaid to illustrate the 
gradual changes in galaxy flattening and orientation. Panels (D) 
and (E): Residual image (observed — model) showing a weak resid- 
ual bar, with a peak intensity of ~ 0.02 e pixel" 1 , corresponding 
to /ig ~ 23.8 mag arcsec -2 . 

spiral arms and non-axisymmetric bending modes — an 
attractive feature of the 2D method. However, full galaxy 
decomposition is, in practice, not always straightforward, 
particularly for nearby galaxies observed at HST/ACS 
resolution. In many situations, it is not clear how many 
components arc needed to fully fit a galaxy, and the phys- 
ical origin of each component may not be obvious. For 
example, sometimes multiple surface brightness profiles 
are required to fit what may be the same photometric 
component (see iPeng et "all [20021 for examples) due to 
the fact that the models used in 2D methods have fixed 
center, ellipticity, and position angle, and have difficulty 
characterizing a galaxy profile in which these parameters 
are not intrinsically constant on all scales. 

The method of ID profile fitting used in this work, 
however, allows the aforementioned parameters to vary, 
and we are therefore usually able to cleanly fit an en- 
tire galaxy with a ID model. A demonstration of this 
is shown in Figure I21[ where we examine the structure 
of FCC 190 (panels A-C), and plot the ellipse model 
derived from our fitted elliptical isophotes (panels D-E) . 
This figure illustrates the striking change in the shape of 
this galaxy's isophotes when moving from small to large 
scales, and how this effect is well-captured by the model. 
The residuals of the fit (panel F) are relatively clean, and 
reveal a weak (fx g ~ 23.8 mag arcsec -2 ) central bar. To 
compare to a 2D fit, the inner 10 x 10" residuals from 
fitting IS and 2S profiles to FCC 190 using GALFIT are 
shown in Figure 125b . Clearly, two Sersic profiles with 
fixed ellipticity and position angle are unable to fully 
parameterize this galaxy. However, the penalty in this 
approach is that the information about the shapes, sizes, 
and relative position angles of various galaxy components 
is lost, as their surface brightness profiles blend together 
into a single component which describes the galaxy on 
global scales. 

Our study is concerned with the properties of the nu- 
clei in comparison to their host galaxies, and with the 
global trends in these properties as a function of galaxy 
luminosity or mass. Thus, we are not interested in a 
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Fig. 22. — Top left: Values for nucleus g-band luminosity obtained from the ID and 2D fits (ordinate and abscissa, respectively). The 
filled blue squares show the host galaxies, the filled black circles indicate the nuclei from galaxies fit well by a S2 profile in 2D, and the 
filled magenta circles represent the nuclei belonging to galaxies in which more than two Sersic components (i.e., S>2) were required in 2D. 
The dotted red line is the relation where the parameters obtained from both methods are equal. Galaxies and nuclei where the measured 
magnitudes differ by more than 0.5 mag between methods are labeled. The root mean squared (RMS) error around the magnitude sample 
mean is shown for the galaxies, nuclei, and nuclei again after applying 3<r-clipping. Top right: Same as the top left, except for Sersic 
indices. Galaxies and nuclei where the measured Sersic indices differ by more than 1.0 between methods are labeled. We note that labels 
for the blue filled squares are to the left of the points, while those for the magenta open squares are to the right. Bottom: Same as top, 
except for galaxy (left) and nucleus (right) geometric mean effective radii. Galaxies and nuclei where the measured effective radii in ID 
and 2D differ by more than 0.1 in the log are labeled. The black dashed line shows the best-fit line to the galaxy effective radii, with a 
fixed slope of unity. 



full decomposition of any large-scale galaxy structure; 
rather, we are seeking to characterize the main galactic 
body component as a whole, so ID techniques are most 
appropriate for our study. However, it is important to 
ensure that the nucleus structural parameters extracted 
using ID methods are robust. To test this assertion, we 
perform surface brightness profile fitting in 2D, and com- 
pare the results obtained using the two techniques. 

A.l. Procedure 

To perform our 2D analysis, we use GALFIT 
(|Peng et all 120021 . l20ll . a program that fits galaxy im- 



ages using multi-component 2D intensity profiles, using 
an iterative downhill gradient Levenberg-Marquardt al- 
gorithm. This 2D analysis is performed on all galaxies in 
our sample with Bt > 13.5, a cutoff which was chosen 
to include most of the nucleated galaxies, while avoid- 
ing those that arc much more challenging to fit in either 
ID or 2D. Galaxies brighter than this are known to of- 
ten show a complex structure, regardless of their clas- 
sification as Es, SOs, dEs or dSOs. For instance, some 
of the brighter dEs are known to show subs tructures in- 
cludin g disk s, spiral arms, and bars (e.g., iLisker et al.l 
l2006bLl2007h . Likewise, more massive galaxies — often 
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Fig. 23. — Top: The difference between 2D and ID galaxy mag- 
nitudes against ID galaxy magnitudes. Bottom: Same as the top, 
but for nuclei magnitudes. 



classified as Es and SOs 



freque ntly show similar mor- 
phological complexities (see, e.g.jBender fc Moellenhofl 

i no^L /~i l Ij. _l 111 r\r\nL Iivt: ' "^j^^T^^onTTcTTTTT^^^^ 



1987t Combes et al. . f990; N ieto et al.lll992HScorza et al 
19981; iFerrarese et alj|2006bt IKrainovic et al.H2011l) . The 
substructures identified in these early-type galaxies could 
either be a sign that they are intrinsically more complex 
objects, or a selection effect arising from their higher 
luminosities and surface brightnesses, which aid in the 
detection of these distinct components. In any case, the 
sample of galaxies used in our 2D analysis consists of 27 
galaxies, 24 of which are found to be nucleated in our 
ID analysis. This sample includes roughly equal num- 
bers of galaxies listed in Table 2 of iPaper J as "giants" 
(E/SO) or "dwarfs" (dSO, dE, etc), although such clas- 
sifications such be viewed with caution since there can 
be significant disc repancies among classifiers: see, e.g., 
IChen et al.l (|2010f ) and Paper III where issues relating to 
the morphology of ACSVCS and ACSFCS galaxies are 
explored in more detail. 

Our analysis proceeded by first measuring the back- 
ground sky value. To do so, we used SExtractor to mask 
out any background sources, and then convolved this 
mask with a Gaussian in order to thoroughly cover any 
diffuse outer edges. The galaxy was then masked with 
an ellipse of geometric radius length between five and six 
effective radii (determined from the ID analysis). We 
then used the biweight value of the remaining pixels as 
the sky value for each of the four ACS chips. Although 
the sky value between different chips was found to vary 
up to ~ 1 count, we found that such a count difference 
resulted in no more than 5% difference in any of the fit- 
ted parameters; we therefore adopted the average of the 
biweight estimates for each of the four chips as the sky 
value. 

We began by fitting each galaxy with a single Sersic 
(SI) profile. We then attempted to fit each of the 24 
galaxies classified as nucleated in ID by adding a second 
Sersic component (for the central nucleus). In 13 cases, 
it was possible to fit the nucleus with a Sersic model 
with all fit parameters varying freely. Five more galaxies 
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Fig. 24. — Panel (a): GALFIT residuals from an SI model fitted 
to FCC 152, showing the inner 10" X 10" region. Panels (b)-(c): 
Same as for panel (a) but for FCC 143 and FCC 26. 

required a prior on the nucleus Sersic index which was 
fixed at n — 2 in analogy with Galactic GCs. For the six 
remaining galaxies, GALFIT was not able to converge on 
a nucleus with only one Sersic component fitted to the 
main body of the galaxy; at least one other large-scale 
component needed to be added for before a fit to the nu- 
cleus could be achieved. However, in one case (FCC 43) 
the nucleus parameters were flagged as having caused nu- 
merical convergence issues, and thus we do not include 
the S>2 fit in our results. In all cases, we did not impose 
any constraints on the relative position angles of of the 
fitted components. 

The above procedures are summarized in Figl22| where 
we have plotted the ID versus 2D magnitudes, Sersic in- 
dices, and effective radii for the galaxies and nuclei from 
our sample. For the galaxies that require more than two 
Sersic components to fit the nucleus, we use the param- 
eters from our 2D SI fit to plot the galaxy portion. Al- 
though the galaxy main body (filled black circles) param- 
eters are generally in good agreement from both tech- 
niques, we note a slight offset in effective radius, where 
those obtained from the 2D fits are usually somewhat 
smaller than in ID (by a factor of 0.94 ± 0.02, derived 
from least-squares fit, with a fixed line slope of 1, to 
the galaxy main body effective radii in the log). The S2 
nuclei (filled blue squares) are also relatively consistent 
between techniques, although with some notable outliers 
that will be discussed below. Finally, the non-S2 nuclei 
(magenta open squares) appear to show the most scatter. 
We note that the scatter in nucleus magnitudes appears 
to be the most significant for the brightest nuclei, prob- 
ably due to the increased difficulty of extracting nucleus 
parameters from structurally complex galaxies that of- 
ten have high central surface brightness. This can be 
seen clearly in in Figure [23l where we have plotted the 
magnitude differences as a function of ID magnitudes. 
We now discuss findings for galaxies in these different 
categories. 

A. 2. Non-Nucleated Galaxies (SI) 

There are three galaxies in our 2D sample where we 
do not find a nucleus in our ID analysis, a result with 
which we find full agreement in 2D. It is interesting to 
examine the residuals of a single- Sersic fit to these objects 
individually to determine why they are not found to be 
nucleated, since it is the lack of a nucleus that is unusual 
for galaxies in our sample. 

The residuals of FCC 152 (Figure 124k ) reveal large 
amounts of dust, but with no nucleus-like object present 
in the central regions. FCC 143 (Figure l2~lb) shows a 
small bar in the residuals, which appears to have a bright 
excess in the center £3 Nevertheless, even with the ad- 

22 Performing a ID fit to the residuals, we find R e ~ CK'02 and 
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Fig. 25.— Panel (a): GALFIT residuals from SI (left) and S2 (right) models fitted to FCC 190, showing the inner 10" X 10" region. 
These results should be compared to the ellipse model (ID) results shown in Figure 1211 Panels (b)-(r): Same as for panel (a) but for 
the galaxy labeled. Galaxies have been ordered by increasing blue magnitude (decreasing luminosity) from the FCC. 



dition of one or two more large-scale Sersic profile com- 
ponents, GALFIT is unable to fit a central nucleus. Fi- 
nally, the low-mass galaxy FCC 26 (Figure [Mh) has two 
bright compact objects located 0'.'95 and l'/37 away from 
the galaxy photocenter. However, it is unclear if either 
of these objects in this actively star-forming, "dE/dlrr 
transition" galaxy can unambiguously called a true "nu- 
cleus" . 

A. 3. Nucleated Galaxies Fit With Double- Sersic 
Profiles (S2) 

Of the 23 nucleated galaxies in our 2D sample, we are 
able to fit the galaxy and nucleus using an S2 profile for 
18 systems. The residuals from the SI and S2 fits to 

g ~ 20.19 mag for the inner and R e ~ 0"22 and g ~ 20.63 mag for 
the outer component. However, due to the complexity of the inner 
structure, we consider these results to be uncertain, and retain the 
non-nucleated classification. 



these galaxies are shown in Figure [25l where the galax- 
ies are displayed in order of increasing blue magnitude 
from the FCC. This figure illustrates how the structural 
complexity of the galaxies seems to increase, and then 
diminish, as their luminosity decreases — reaching an 
apparent maximum in the range —19 < Mb < —17 
or 10 10 - 4 > M*/M Q > 10 9 - 6 — with the residuals 
for the faintest galax i es appearing much cleaner (see 
iFerrarese et"afl 120061* ILisker et alJ I2006H ) . Of course, 
part of this apparent simplicity is likely related to the 
lower S/N of the available imaging for the faintest and 
lowest surface brightness systems. 

For five of the S2 galaxies, the Sersic index of the nu- 
cleus needed to be held fixed during the fit. FCC 190 
(Figure^), FCC 55 (Figure[2SU), FCC 95 (Figured!:), 
and FCC 136 (Figure l2"5k ). all have substructure such 
as bars that overlap with the nuclear region which the 
second Sersic component attempts to fit. By fixing the 
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TABLE Al 

ID and 2D Nucleus Parameters for Multi-Component 
(S>2) Galaxies 
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Sersic index of the nuclei at n = 2 (appropriate for Galac- 
tic GCs), GALFIT is able to fit the nucleus, with the 
resulting magnitude and effective radius of the nucleus is 
in agreement with the ID results in all cases except for 
FCC 190, which is discussed below. The other galaxy 
that requires a fixed nucleus Sersic index, FCC 335, con- 
tains a large amount of dust in the central regions, and if 
the nucleus Sersic index is allowed to vary, then the nu- 
cleus effective radius and Sersic index converge to very 
small values that GALFIT warns may cause numerical 
convergance issues which cause the final solution to be 
unreliable. The differences between the ID and 2D re- 
sults for this galaxy are also discussed below, where we 
describe nuclei that are notable outliers in Fig[22] Specif- 
ically, the nuclei of FCC 190, FCC 335 and FCC 90 have 
ID and 2D magnitude differences of > 0.4 mag and frac- 
tional differences in their effective radii of > 0.5. 

FCC 190: This nucleus is 0.55 mag brighter and twice 
as large in effective radius in the 2D fit than in ID. 
The galaxy exhibits distinct "peanut-shaped" residuals, 
as shown in Figure [25k . It should be noted that after 
fitting a both a bulge and a disk component along with 
the nucleus, the nucleus magnitude and radius are still 
notably disparate. 

FCC 335: In 2D, the nucleus is 0.42 mag brighter, but 
four times smaller, than in the ID fit. The 2D residuals 
are shown in Figure [25b . This galaxy has a large amount 
of dust, and the center was held fixed during the ellipse 
fitting for the ID analysis. However, the position of the 
ID fit ellipse center is actually ~ 0.5 pixels away from 
the nucleus (as determined by GALFIT and confirmed by 
eye). This could cause the ID analysis to overestimate 
the nucleus effective radius and underestimate the mag- 
nitude, as the light from the nucleus effectively becomes 
smeared out. 

FCC 90: This nucleus is 1.12 mag brighter in 2D than 
in ID. The residuals of a single-Sersic fit (Figure 
show a bright central nucleus as well as a secondary 
fainter object ~ C/25 away. This second object is the 
cause of a small secondary bump in the ID surface bright- 
ness profile (see Figure [3]). After simultaneously fitting 
this secondary object, the nucleus is still found to be 
1.08 mag brighter in 2D than in ID. Like FCC 90, there 
is large amounts of dust in the center of this galaxy, and 
the center was held fixed for the ID ellipse fitting, at a 
point ~ 1 pixels away from the 2D nucleus center, which 
may partly account for the smaller and brighter nucleus 
found in 2D. 

A. 4. Nucleated Galaxies with Multiple Large-Scale 
Components (S>2) 

There are six nucleated galaxies in our sample which 
we were unable to model in 2D using an S2 model, as a 



second Sersic component in GALFIT will, even with its 
Sersic index held fixed, attempt to fit a different com- 
ponent of the underlying galaxy. Thus, we need to add 
a second, or even third, Sersic component to the main 
body of the galaxy in order to fit the nucleus (i.e., three 
or four components in total). A comparison of the ID 
and 2D nucleus parameters for these six galaxies is given 
in Table IA1I Note that the nuclei in several of these 
galaxies are candidates for "hybrid nuclei" with a com- 
plex structure (e.g., a compact, high surface brightness 
component embedded in an small-scale disk-like feature) 
that are suggestive of multiple, parallel formation pro- 
cesses (see fllXTl and 14.3.21) . 

FCC 43: A small (< 10" in diameter) round central 
component, seen in Figure [25k needs to be fitted before 
it is possible to fit to the nucleus. However, because 
of the small size of the fitted nucleus (i.e., an effective 
semi-major axis of 0.39 pixels), the output of GALFIT 
indicated that this parameter may have caused numeri- 
cal convergence issues, making all parameters from this 
solution unreliable. We therefore did not include these 
results in Fig |2"21 

FCC 310: After fitting with a single Sersic profile, the 
bar and and outer envelope-like structure of this galaxy 
become apparent from the residuals, seen in Figure [2"5"b . 
To fit the nucleus, we must first fit an n — 2.38 bulge- 
like component, an n — 0.26 bar, and an n = 0.20 outer 
envelope. After these three components are fit, the nu- 
cleus appears quite clearly in the residuals, and it can 
be fitted by adding a fourth Sersic profile with its Sersic 
index held fixed at n = 2. 

FCC 249: A single-Sersic fit reveals a peanut-shaped 
residual in the center (Figure [25b ) , with a possible nu- 
cleus. After a second small component is added, a nu- 
cleus becomes apparent in the residuals. The nucleus can 
then be fitted with its Sersic index held fixed at n = 2. 

FCC 148: This galaxy shows a very boxy inner bulge, 
with X-shaped isophotes in intermediate regions (Fig- 
ure [25b0 . Since we are unable to fit the host and nucleus 
with a double-Sersic model, a second large-scale compo- 
nent with disk-like properties (n = 1.04, and an axis ratio 
of 0.36) is added, after which GALFIT will converge on 
the nucleus. Although a nucleus is not very prominent 
in the two-component fit residual, the Sersic index of the 
bulge- like component grows to n = 9.35 if a nucleus is 
not included in the fit. After a nucleus is included, the 
fitted bulge Sersic index is n = 5.1. The disk-like com- 
ponent does not change significantly with the addition of 
the nucleus. 

FCC 301 : The complex structure of this galaxy, seen 
in the single Sersic component fit residuals in Figure [25b . 
can be appreciated from the ID surface brightness profile 
(Figure[3]), where the intensity is slightly over-subtracted 
at 1", and then under-subtracted out to ~ 5". There 
are also bright outer wings, at > 10" scales. After a 
single-Sersic component is fitted, a second component 
will converge on the larger bright central excess, and a 
third Sersic profile will then fit the nucleus. However, 
the bright residuals show that the the main body of this 
galaxy is not well described in 2D, even by two Sersic 
profiles. 

FCC 204 : As is the case for FCC 43, there appears to 
be an embedded disk in this galaxy, which can be seen 
in the residuals of a single Sersic fit, Figure l25"f . The 
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(a) 




FCC 204 m m « 

Fig. 26. — -Row ('aj: The first three panels show, from left to right, 80" X 80" regions centered on FCC 43 with GALFIT residuals from 
an SI fit, a two-Sersic component fit, and a two-Sersic component plus Sersic nucleus component fit. The last three panels show the same, 
but magnified to show the inner 10" X 10". Rows (b)-(f): The same as panel (a) but for the galaxies labeled. In the case of FCC 310, the 
residuals show the results found using three, rather than two, Sersic components. 



nucleus is found to be slightly fainter in the 2D fit than 
in ID. It is possible that in ID, the central disk might 
be contributing to nucleus component and causing the 
nucleus luminosity to be over-estimated. 

Overall, the 2D nuclei parameters from these complex 
galaxies are in reasonable agreement with those found in 
the ID analysis. In terms of magnitude, only FCC 148 
and FCC 249 show differences of > 0.5 mag, although all 
except for FCC 310 show discrepancies of > 50% in effec- 
tive radius and Sersic index. However, the differences do 
not appear to be systematic, in the sense that there does 
not seem to be consistent under- or over-estimation of a 
specific parameter in ID or 2D. The nuclei parameters for 
these cases are likely to be more uncertain overall, and 
thus larger differences between the extracted parameters 
in these structurally complex galaxies is to be expected. 



A. 5. Nuclei Ellipticities 

One advantage of performing the 2D analysis is that 
we are able to measure the nuclei ellipticities. In general, 
we find the nuclei to be flattened, with median elliptic- 
ities of 0.20 and 0.25 for the full and only S2 samples, 
respectively. The top right panel of Figure [27] shows a 
histogram of the ellipticity distribution, while in the top 
left panel, we plot ellipticity against host galaxy magni- 
tude. A least-squares fit to all of the nuclei hints at a 
trend of increasing nuclei ellipticities with galaxy lumi- 
nosity, 

£nuc = -(0.062 ± 0.027) B T + (1.2 ± 0.40), (Al) 

and after removing the more uncertain S>2 nuclei, we 
obtain a significant relation: 

e„uc,S2 = -(0.10 ± 0.03) B T + (1.8 ± 0.45). (A2) 
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Fig. 27. — Top left: Nuclei ellipticity vs. host galaxy magnitude, 
for the S2 (closed black circles) and S>2 (open magenta circles) 
nuclei. The dotted and dashed lines show fits to the full sample 
and only the S2 nuclei, respectively. Top right: Histogram of the 
nuclei ellipticities. The gray regions represents the S2 nuclei, while 
the magenta regions indicate S>2. Bottom: Nuclei ellipticity plot- 
ted against the absolute difference between nuclei and host galaxy 
position angles. 

This result suggests that nuclei in brighter (and higher 
mass) galaxies are more flattened, and may be more likely 
to contain edge-on disk-like components. 

In the bottom panel of Figure [27J we show nuclei el- 
lipticities versus the difference in position angle between 



the nuclei and their host galaxies. We find that the most 
highly flattened nuclei are aligned with their hosts, al- 
though over one third (7/18) of our sample are signifi- 
cantly (AP.A. > 20 deg) misaligned. 

A. 6. 2D Analysis Conclusion 

Figure [55] shows that there is reasonable agreement 
between the nuclei parameters measured in ID and 2D. 
We conclude that the ID nucleus parameters are for the 
most part robust, and note that the brightest and most 
structurally complex galaxies — which typically have 
Ms(l") 19 ma § arcsec -2 — present a challenge for 
measuring nuclei parameters using either approach. 

Indeed, even in cases where adding a second or third 
profile to the main body is required to fit the nucleus in 
2D, it is unclear how many components must be added 
until a "best" fit is actually achieved, and it is usually 
difficult to say whether one method yields parameters 
closer to those of the true nucleus. In our study of For- 
nax nuclei, we are primarily interested in extracting the 
nuclei parameters relative to the average outer profile. 
Although much of the power of 2D techniques lies in 
their ability to fit multiple large-scale components, in 
galaxies that require more than one outer Sersic profile, 
it becomes more difficult to perform a fully objective and 
homogeneous comparison between the nuclei and galaxy 
parameters. Our analysis therefore uses the results from 
our ID fits, which meets the dual criteria of objectivity 
and homogeneity. The general consensus between meth- 
ods indicates that the main conclusions in this work are 
independent of the specific approach used to model the 
galaxies and their nuclei. 
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